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Abstract 

We investigate the relationship bet'ween colour confinement and the gauge independent Cho-Duan- 
Ge Abelian decomposition. The decomposition is defined in terms of a colour field n; the principle 
novelty of our study is that we have used a unique definition of this field in terms of the eigenvectors 
of the Wilson Loop. This allows us to establish an equivalence between the path ordered integral 
of the non-Abelian gauge fields with an integral over an Abelian restricted gauge field which is 
tractable both theoretically and numerically in lattice QCD. We circumvent path ordering without 
requiring an additional path integral. By using Stokes’ theorem, we can compute the Wilson 
Loop in terms of a surface integral over a restricted field strength, and show that the restricted 
field strength may be dominated by certain structures, which occur when one of the quantities 
parametrising the colour field n winds itself around a non-analyticity in the colour field. If they 
exist, these structures will lead to a area law scaling for the Wilson Loop and provide a mechanism 
for quark confinement. Unlike most studies of confinement using the Abelian decomposition, we 
do not rely on a dual-Meissner effect to create the inter-quark potential. 

We search for these structures in quenched lattice QCD. We perform the Abelian decomposi¬ 
tion, and compare the electric and magnetic fields with the patterns expected theoretically. We 
find that the restricted field strength is dominated by objects which may be peaks a single lattice 
spacing in size or extended string-like lines of electromagnetic flux. The objects are not isolated 
monopoles, as they generate electric fields in addition to magnetic fields, and the fields are not 
spherically symmetric, but may be either caused by a monopole/anti-monopole condensate, some 
other types of topological objects or a combination of these. Removing these peaks removes the 
area law scaling of the string tension, suggesting that they are responsible for confinement. 

Keywords: Quantum chromodynamics. Lattice gauge theory. Confinement of Quarks 
PACS: 12.38.-t, 12.38.Aw, 11.15.Ha 


1. Introduction 


An enduring problem in QCD is to find the mechanism which causes quark confinement, which 
is known to be non-perturbative in its origin. Although several models have been pr^osed - for 
example, center vortices [l|, and a dual Meissner effect due to magnetic monopoles i i i i - 
there has not yet been a convincing demonstration that any of them are correct. In this initial 
study, we investigate the Cho-Duan-Ge (CDG) Abelian decomposition (sometimes referred to as 
the Cho-Faddeev-Niemi decomposition) i 0, i i [13, and concentrate on one of the signals 
of confinement, a linear static quark potential. Unlike Dirac and ’t Hooft (Maximum Abelian 
Gauge) monopoles, the CDG decomposition allows for monopole solutions while respecting the 
gauge symmetry and does not require a singular gauge field or an additional Higgs field. The 
decomposition is constructed from a colour field, n, which may be built from a SU(fVc) matrix 
9, where S\J{Nc) is the gauge group of QCD. Each n dehnes a different decomposition, and 
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an important question is the best way of choosing this field. Recent work [ll|, [l^ [l^ [l^ has 


demonstrated that the magnetic part of the field strength dominates the confining string, using a 
decomposition constructed from one particular choice of 0 G SU{Nc)/U{Nc — 1); however in this 
case only one of the possible Nc — 1 types of monopole is visible. 

Here we consider a different choice of 0 G SU{Nc) /. Our initial goal was to 
investigate whether monopoles apparent in this construction may also lead to confinement^. Con¬ 
centrating on the Wilson Loop, an observable used to measure the string tension, we show that 
the path ordering may be removed by diagonalising the gauge links along the Wilson Loop by a 
SU{Nc)/{U{l))^<^-^ field 9. This, in principle, allows us to use Stokes’ theorem to express the 
Wilson Loop in terms of a surface integral over an Abelian restricted gauge field strength tensor. 
We use the CDG decomposition, constructed from a colour field for a diagonal Gell- 

Mann matrix , to find a consistent construction of this restricted field across space-time, and not 
just for those gauge links along the Wilson loop where it was originally defined. Our relationship 
for the string tension in terms of this restricted field is exact: we do not require any approximations 
or additional path integrals. We search for topological structures in this field strength, and find 
that discontinuities in the colour field related to the winding of one of the parameters describing it 
around a particular location may lead to large ‘electric’ and ‘magnetic’ fields (deriving the terms 
from an obvious analogy to electromagnetism) in both the restricted field strength and the part of 
the field strength derived from the colour field; in particular the component of the field strength in 
the plane of the Wilson Loop may exhibit structures, decaying with a power of the distance in all 
directions, which would lead to an area law scaling of the Wilson Loop. We label these structures 
merinthons. We also examine the other components of the electromagnetic field strength and find 
that if certain simplifying assumptions are reliable they will be dominated either by similar peaks 
or lines of electromagnetic flux. We do not show that this winding occurs in practice, but if it 
does it provides an explanation for the confining potential. We also discuss how string breaking 
might arise in this picture at large distances. 

The second part of this article checks whether these structures are present using SU (3) quenched 
lattice QGD. We have generated 16^32 and 20^40 gauge field configurations, and used them to 
test this theory. We have concentrated on gauge invariant observables to avoid certain ambiguities 
with the gauge fixing: the restricted field strength is gauge invariant, but our choice of 9 field is 
not, making it harder to directly observe merinthons in the underlying field. We perform the 
Abelian decomposition and calculate 9, n, the field strength tensor and the portion of the field 
strength tensor due to the structures in the 9 field. Our goal in this initial study is to show that the 
Electromagnetic field contains the expected structures, and that the merinthons can account for 
the string tension. Our simulation is limited because although in principle we should recalculate 
9 and use a different Abelian decomposition for each Wilson Loop, we do not do so due to limited 
computer resources; nor have we considered the deconfinement transition, different representations 
of the gauge group, flux tubes of the confining string, or the relationship between this picture of 
confinement and the center vortex or dual Meissner models. These would have to be considered 
in later works. 

. In section [5] we discuss the Abelian 


Preliminary results were presented in |16l |. 17| and 


decomposition and its relation to the Wilson Loop and thus quark potential. We discuss the 
differentiability of the colour field in section [3l and then how confinement may arise in this picture 
in section m After a brief note on string breaking in section [SJ we present numerical evidence in 
section [6] and conclude in section [71 There are four appendices outlining the details of some of the 
calculations and our numerical methods. 


different choice of a SU{Nc)/{U{l))^'^ ^ Abelian decomposition was described in m , but without a 
discussion of the relationship to quark confinement. 
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2. The Abelian Decomposition and Stokes’ theorem 

The Wilson Loop in an SU(Ac) gauge theory is defined as 

i\c 

for a closed curve Cs of length L which starts and finishes at a position s, where V represents 
path ordering and the gauge field, can be written in terms of the Gell-Mann matrices, A“, as 
2 A“ A“. We will use the summation convention that the superscripts a, 6,... on a Gell-Mann matrix 

implies that it should be summed over all values of a, A“A“ = X]a=i A“A“, while the indices 
j, k,... are restricted only to the diagonal Gell-Mann matrices, so, in the standard representation, 
A-^A-’ = J2j=3,8,...,N^- _^A'^ AL We shall often leave the gauge field dependence of W and Wl 
implicit. 

The Wilson Loop when Cg is an ii x T rectangle, with spatial extent R and temporal extent 
T, can be used to measure the conhning static quark potential, V{R), 0 

T(i?) = - hm log((W4a]))/r, (2) 

1 —¥oo 

where (...) denotes the vacuum expectation value. It is expected, and observed in lattice simu¬ 
lations, that for intermediate distances the confining potential is linear in ii, so V (R) ^ pR -I- c, 
where p is the string tension, and c is a constant. At very small distances, the potential is ex¬ 
pected to be Goulomb, while in the presence of fermion loops at very large distances the string is 
broken and the potential becomes independent of R [i^. The main focus of this work is on the 
intermediate regime, so we expect that the expectation value of the Wilson Loop will scale with 
the spatial and temporal extents of the Wilson Loop as 

{Wl[Cs\)-- e-P^^. ( 3 ) 

This is known as the area law scaling of the Wilson Loop. As discussed in [Appendix B[ this is 
satisfied, if on each individual configuration, the Wilson Loop scales as 

(4) 

where F is randomly distributed from configuration to configuration (according to one of a certain 
set of distributions which includes those discussed in this work) with a mean value proportional 
to the area contained within the curve Cs- The eventual goal (and this work is intended as a step 
towards that goal) is to demonstrate from first principles that equation (|4]) is satisfied in pure 
gauge QCD (and later full QCD), and thus that the quarks are linearly confined. 

A difficulty with evaluating equation o is the path ordering. If the fields Afj,{x) at different 
X and p commuted with each other, then we could ignore the path ordering, and use Stokes’ 
theorem to convert the line integral to a surface integral. The problem would then reduce to 
showing that there was some flux flowing through the surface so that the surface integral was 
proportional to the area (perhaps by counting lines of flux). However, A^ are non-Abelian fields: 
we cannot immediately do this. Our approach is then first to construct an Abelian field A^(a;) 
so that Wl [Cs , A] = Wl [Cs , A], and the calculated string tension does not depend on which of 
the two fields we use. We may then remove the path ordering, apply Stokes’ theorem to replace 
the line integral with a surface integral over some field strength, and then show that the surface 
integral is proportional to the area enclosed within the loop. 

First we shall define what is meant by path ordering. We split Cs into infinitesimal segments of 
length 5(7, and define the gauge link as U„ G SU{Nc) = V[e-^P .Q<a<L 

represents the position along the curve and we write A^ = A^(ct) (a:(cr)). We have assumed and 
will require throughout this work that the gauge field, A, is differentiable. This limits us to only a 
certain subset of gauges, and once we have found a suitable gauge we are restricted to continuous 
gauge transformations, i.e. only those gauge transformations which can be built up from repeatedly 
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applying infinitesimal gauge transformations, —>■ + g~"^d^a + i[a,A^], where a = a“A“ and 
dfj,a are both infinitesimal. We also neglect the effects of the corners of the Wilson Loop; the 
discontinuity in Aa- at the corner can, for example, be avoided by using a rounded corner. The 
path ordered integral over gauge fields is defined as the ordered product of the gauge links around 
the curve in the limit Jcr —> 0. 

W[Cs] can be written in this lattice representation as 

L — 5a 

W[C.]=lim n (5) 

5(7^0 ^ ^ 

cr=0,(5cr,2(5(7,... 


We wish to now replace Ug by an equivalent Abelian field. 

Previous work (^.[^ , [ll | has achieved this by inserting the identity operator, as parametrised 
below, between each pair of gauge links, 


I DP 


d9 


( 6 ) 


where 9 G SU{N), e^c is a unit vector in colour space, and d9 is a suitable Haar measure for 


9. 


Proof that this is an identity operator is given in [ll|. This allows one to replace Ug with 
Ug = +r9dae )ejvc _|_ o{5a‘^). The exponent is Abelian, so there is no need for path 

ordering. In effect, this replaces the path ordering with a path integral over 9. It can be shown that 
the exponent is related to the CDG Abelian decomposition, which allows for certain monopole 
solutions which reside in the field 0, and these monopoles may lead to confinement via a dual 
Meissner effect. 

We believe that this approach contains a number of difficulties which have to be resolved before 
it can be used practically: 

1. It introduces an extra gauge degree of freedom to the system - contained within the variable 

2. It is necessary to evaluate an additional path integral, making calculations harder; 

3. In gauge groups larger than SU(2), only a small subset of the possible CDG monopoles are 

visible to the restricted field (i.e. only the monopoles contained within the colour 

direction, n will be defined in equation (O); 

4. There is no obvious relation between the physical gauge field and the unphysical 9 field which 
contains the monopoles and leads to confinement. We expect topological features within the 
gauge field to contribute to confinement; and this is only possible if the monopoles in the 
9 field are related to the monopoles in the gauge field. However, given that we integrate 
over every possible 9, it is difficult to see how the gauge field itself can affect the confining 
potential if both the CDG monopole picture and this model are correct; 

5. It is possible to introduce a gauge transformation which leaves 9 constant across space time, 
destroying any monopole structure; 

6. The 9 field at each location along the curve is independent (as each needs to be integrated 
separately), and therefore e^j^^{9Ag9^ + i9dg9^)eNc will not be a smooth function of the 
position for a general configuration of 9. This makes the conversion lim^CT-j-o X] I 
problematic - this is required in the conversion from a product of gauge links to the exponen¬ 
tial of a line integral of Ag around the Wilson Loop - and there are also difficulties defining 
dg9. Some regularisation procedure would be required to ensure that 9 is smooth. Normally, 
in a path integral in Euclidean space time, discontinuities if the fields are suppressed by the 
differential terms in the action, so we can rightly assume that the fields are differentiable. 
However, in this case 0_^es not contribute to the action, so there is no similar suppression 
of discontinuities (see 


23| for a more detailed criticism along these lines). 


In [ll|, many of these challenges are resolved by fixing the additional degrees of freedom by a 
procedure analogous to gauge fixing in QCD. This, in effect, selects one particular 9 field and 
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uses that for the subsequent analysis (an alternative way of expressing the approach of m is to 
say that their procedure fixes various components of the gauge field but this breaks gauge 
invariance). However, this breaks the identity in equation ([S]) and thus the relationship between 
the Wilson Loop of the restricted gauge field U and the Wilson Loop of the original gauge field 
U. In practice, to pursue a program along these lines would require sampling numerous different 
choices of 0-fixing. 

In view of these challenges, we propose that different approach should be taken when trying 
to express the Wilson Loop in terms of an Abelian field. 

We introduce a field 6^ = 0{x{a)), which, for the moment, we shall take to be an element of 
U(fVc), at each location along Cg and insert the identity operator 0^^^ between each of the gauge 
links. 0 is chosen so that O^^Uada+Sa is diagonal. There are L/Sa gauge links along the path, and 
we introduce L/6a 0 fields, so there is no obvious reason why the system cannot be solved. In 
fact, it is easy to construct a solution: it is easy to show (see [Appendix AT|) that 0s contains the 


eigenvectors of IT[Cs]: W[C's]0s = 0se*^^^ diagonal “ ^ for some real . 

As the phases of the eigenvectors are arbitrary, this definition only determines 0 up to a 
([/(l))^*^ transformation 0 —^ 0%. x makes no difference to any physical observable, but for prac¬ 
tical purposes it is useful to select the phases and ordering of the eigenvectors by some arbitrary fix¬ 
ing condition to give a unique choice of 0 S SU{Nc)/{U. We define SU{Nc)/ 
by considering the following parametrisation of a [/ {N) matrix: 


/ cos tti 

i sinoie*'’^ 

0 . 

■'1 


( cos 02 

0 

f sin 026*'’^ 


i sinaie“*^i 

COSOi 

0 . 


0 

1 

0 

0 

0 

1 . 



i sin 026“*'’^ 

0 

cos 02 


V 



■■ ) 


1 : 



) 




XJ diagonal 


dj ) 


tti, Ci, c?o and dj are real parameters, and there are Nc{Nc — l)/2 Givens matrices (i.e. one for 
each of the possible ways of embedding a 2 x 2 matrix into a Nc x Nq matrix) parametrised by 
one particular a and c. An SU{Nc )/matrix is parametrised in the same way, but 
without the final (17(1))^'^ term (i.e. by setting dj and do to some arbitrary fixed value, most 
conveniently dj = do = 0). 

Under a gauge transformation, —>■ ^crUcrJ^„+&a ^ G SU{Nc)^ 0 —t A0x, 

where the ([/(l))^"^"^ factor x depends on the fixing condition (in this case A0 G SU{N) so 
there is no contribution to x from do)- This follows from the definition of 0 as containing the 
eigenvectors of the Wilson Loop. The Wilson Loop transforms under a gauge transformation 
as W[Cs,U] ^ AgIU[Cg,C/]At, so the operator which diagonalises it transforms according to 
0s —t As0s; although we also need to reselect the 17(1)^'^“^ factor so that the fixing condition 
remains satisfied. With 9\Ua0a+S(y = diagonalfor real u, 


0jlU[C's]0s = fcs ^ 


(7) 


removing the non-Abelian structure and the path ordering without introducing an additional path 
integral. 

Our goal is to apply Stokes’ theorem to convert this line integral into a surface integral, and 
this requires extending the definition of 0 and across the surface bounded by Cg- In practice, we 
construct these fields across all of space time. To generalise 0, we construct nested curves, C*, in 
the same plane as Cg and then stack these rectangles on top of each other in the other dimensions, 
so that every location in Euclidean space-time is contained within one and only one curve. We 
then define 0 so it diagonalises the gauge Links (and only the gauge links) which contribute to 
W[C\U]. 

We cannot naively extend across all of space-time, because its definition requires that all 
the gauge links U are diagonalised by 0, not just those that contribute to the Wilson loop, and in 
general this cannot be satisfied. Instead, we replace the gauge links U with a field U, defined in a 
consistent way so that it is both diagonalised by 0 across all of space time, and equal to U along 
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the path Cg- The first of these conditions means that 

[A^ =0, (8) 

for each diagonal , which can be re-written in the form 

- < =0 < (9) 

This condition is satisfied across all of space-time and for all directions fji. Note that n-1 is inde¬ 
pendent of the choice of x- As we shall see later, the 0-dependence of the restricted field strength 
Ffii, [A] only appears within , and objects contributing to the restricted field strength drive con¬ 
finement. This is the justification of our earlier statement that the choice of x does not affect the 
physical observable, which is the restricted field strength. To give this U field a physical meaning 
we need to relate it to the gauge field U, and we do so via a second field X defined according to 

= U^{x)Ul. (10) 

For later convenience (equation (1501) 1. we restrict X^ by imposing the condition 

trK(Xt,-l^,,)]=0. (11) 

Under a gauge transformation n transforms as nx —>■ AxUxXl (which follows from the transfor¬ 
mation rule for 9) and the requirement that equations (jO]) and dill) are satisfied in every gauge 
leads to the transformation rules Ufj.{x) —>■ AxU^i x^x+iiSa a; AxX^xA\.. Equations 

® and CD are the lattice versions of the defining equations of the gauge independent CDG 
decomposition i0,iS[i3 , which in the continuum is described bjH 



(12) 

D^[A]n^ =0 

(13) 

o 

II 

(14) 

Dfj_[A]a =d^a - ig[A^, a] 

(15) 


(16) 
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with 

(17) 

Equation (HU is the naive continuum limit of equation (US; equations m and (113 are the 
naive continuum limit of equation ([^, and equation m is the continuum limit of equation (HD- 
Proof that elation (USD solves the continuum decomposition has previously been given in (for 
example) |6|, [tI 1ii| , and is repeated in [Appendix A| 


^The original discoverers of this decomposition prefer to call it a gauge independent decomposition rather than 
gauge invariant decomposition. In these earlier models, the choice of 9 was left arbitrary. The procedure was to 
fix to some arbitrary gauge, and then select some 9. Since the decomposition only depends on the choice of and 
proceeds regardless of which gauge was originally chosen, the decomposition was referred to as gauge independent 
to distinguish it from other approaches which required fixing to particular gauges, such as the decomposition based 
on the Maximum Abelian Gauge. Our work differs in philosophy from the original approach because we do not 
need to perform any gauge fixing to extract our observables, since our key quantities (such as the field strength 
and Wilson Loop) are gauge covariant and thus the corresponding observables are gauge invariant. However, our 
particular choice of 9 is gauge dependent, and thus the quantities which we use to parametrise it can only be 
examined after gauge fixing to some arbitrary gauge. The authors of gauge fixed to an arbitrary 

gauge and then performed an Abelian decomposition; we decompose and then (if required, which isn’t the case for 
any physical observable) gauge fix. 
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The condition m may then be interpreted as requiring that tr(n-^ A) = tr(n^ A), the component 
of the gauge held A parallel to n is fully contained within A. In the continuum, the solution for A 
and X is unique, but on the lattice we found that there were sometimes several distinct solutions 
to equations m and (HU. In this case, we choose the solution which had the largest value of 
tr(X), a condition which is both gauge invariant and satished along Cs where U = U and thus 
1 = 1 . 

The continuum dehning equations ensure that the held strength associated with A, dehned 

by 

[D^[A],D,[A]]a = -ig[FM,o\ ( 18 ) 

for any held a in the adjoint representation of the gauge group, satishes F^v[A\ = for some 

real scalars The proof of this follows by substituting each of the helds in turn in place of 
a, using equation (HI to show that [F,n^] = 0 for all the , and noting that the only objects 
which commute with each of the are the other helds. 

We express the restricted held as for real m, and since U = U along 

the curve Cs, we see that IT[Cs, U] = W[Cs, U] = 0sW[Cs,d^U0]0\ = 0se^^ Applying 

Stokes’ theorem to the Abelian held 0l.Ufi^x0x+fLSa gives, if u is differentiable, 

0lW[Cs]0s , ( 19 ) 

=d^ni - (20) 


where F^ (like u) is gauge invariant, S the (planar) surface bound by the curve Cg, and dS an 
element of area on that surface. Note that u does depend on the hxing condition, although F is 
independent of it. 

Equation (j20ll can be easily proved by considering for smooth 0 and u 


Ua = e 




which implies that 

gA^ = -ulnd + i0da0\ 


( 21 ) 


( 22 ) 


[A] —g ^dfj^Aj^ dyA^ t[A^, A^/]^ 

= — uidf^nd — n^df^ul + — uj^[nd ,0dy0^] + ui[n^, 0dfi0^] 

- i0df,0^0dy0'< + i0dy0^0df,0^ + i[0df,0\0dy0^] 

=nddyuj^ - ndd^ul, (23) 

with F^y[A] = ^tT{ndFfiy[A]). We have used the result 

dfj,n^ = [nd ,0df,0^], (24) 


which follows from the dehnition 

We can now consider the Wilson Line around an inhnitesimal plaquette p, which for a smooth 
u held gives 


which leads to 


~ ^x,fiUx+5cril,vUxj^s„y Jjx,v “6 i 


lT[p,6»1'C/6»] = e 


_ iX^{d^ul, — duU^ 


(25) 


(26) 
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and building the integral over the surface bounded by Cg from the product of integrals over these 
small plaquettes, using that the exponent is Abelian, gives equations (ITOl) and (1^ . 

Proof that the definitions of F given in equations dT51) and (1^51) are equivalent is given in 
[Appendix A[ along with a proof that the definition of in equation (1221) is equivalent to that 
given in equation (HU). 

Finally, as alluded to earlier, we note that t)e written in the forms 

F/xi/[A] =^nAdfj.t^n^Ai, - d^tin^A^ - ,0d^9^) (27) 

=T:nAdAT:n^ Ai, - d^tm^ A^ + -^nArinAdfj.n'" ,d^nA), (28) 

2 85 

as is proved in [Appendix A[ equations (IA.25|) and (IA.27I) . These functions solely depend on n, 
and 0 only indirectly through n, and since n is independent of x F is the physical observable 
we want to study, we conclude that the choice of x will not affect any of the physical observables 
we need. 

Equation (HU is only valid if u is differentiable. Equation (|T^ is also similar to what we see 
in QED, which is, of course, not confining. In analogy to QED, we may expect the contribution 
of those portions of space time where u is continuous to have little contribution to the string 
tension. However, we must also add to this equation the effects of discontinuities in u. We do so 
by only extending the area integral over those areas where u is continuous, and add additional 
line integrals around the areas where it is discontinuous. The linear string tension will, at least in 
part, arise from these discontinuities. Since vP is built from the gauge field, and we are working 
on a gauge where this is assumed to be differentiable, and 9, we must therefore consider whether 
9 is differentiable. 

3. Differentiation of eigenvectors 

3.1. SU(2) 

We wish to investigate how 9 varies over space, which requires evaluating the derivative of 
the eigenvectors of the Wilson Loop operator. Suppose that an operator M has two eigenvectors 
(which is the case for SU(2)), tpi and ■02, with eigenvalues Ai and A 2 and after a small change in 

M ^ M + 6M, the eigenvectors change to V'l and ipA with eigenvalues Ai + <5Ai and A 2 + i5A2. 

For normal M (so that the eigenvectors are orthogonal, left and right eigenvectors the same, and 
5M is (anti-)Hermitian or an anti-Hermitian matrix multiplied by M if M is unitary), we may 
parametrise the new eigenvectors as 

="01 cos(/3) + ^/’ 2 e*’^ sin/3 

V '2 ='02 cos(/3) — sin /3, (29) 

and the eigenvalue equations read 

Mipi =Ai'0i 

M'02 =A2'02 

(M + SM){'tfji cos(/3) + V' 2 e*^ sin/3) =(Ai + <5Ai)(V'i cos(,5) + V' 2 e*’’' sin/3) 

(M + SM){'tp 2 cos(/3) - V'le*'’' sin/3) =(A 2 + SX 2 ){'tp 2 cos{P) - V'le*’’' sin/3). 


(30) 










which, after applying and ipl to these equations gives 

5Ai ='ijj\SM'>pi + tan/3 

(5A2 =V'|<3Mi/)2 — iplSM tan ,5 

\/ijj2SMiljiip\SMip2 

tan 2/3 =2-- 

Ai — A 2 + — tplSM'tp2 


At |Ai — A 2 I » 
theory [ 2 ^ 




'ip\dM'4}2 


( 31 ) 


^|(3M'02|, this agrees with the result from first order perturbation 


St/ji 


V'^ - V'i = (1 - 


1 

M -X 


6Mipi. 


(32) 


In this case, the change in the eigenvector is proportional to 6M. However, if the eigenvalues are 
near-degenerate (i.e. the difference between them is of 0(||5M||)), then /3 may be large even as 
i5M —>■ 0, and the eigenvectors would evolve discontinueusly. 

As we will be interested in the spatial derivatives of the eigenvectors of the closed Wilson 
Loop, I will set M to be the Wilson Loop operator W[C'J], and consider the derivatives in two 
directions, one parallel to the Wilson Line, and one perpendicular to it in the plane of the Wilson 
Loop. Any small shift in location in an arbitrary direction can be constructed in terms of a shift 
parallel to the Wilson Loop followed by a shift perpendicular to the Wilson Loop, r is an index 
denoting the size of the Wilson Loop; an increase in r means that we move to the next of the 
nested Wilson Loops used to define 0, so both the spatial and temporal extents of the curve CJ 
would increase by a small amount 2i5r. For an R x T planar Wilson Loop with T > R we have 
r = R/2. We investigate the change of the Wilson loop both when we have a small shift along the 
curve, s —>■ s -|- da, and a small shift perpendicular to the curve, r ^ r + Sr. 

Firstly, there will clearly be discontinuities in 9 when the gauge field is discontinuous, but as 
we require a continuous gauge field for the decomposition to be valid we neglect this case. We 
wish to consider if there can be discontinuities in the eigenvectors of the Wilson line even for a 
continuous gauge field. 

In the direction parallel to the line, we have 

W[c:+s.] = U!W[C:]U,, (33) 

and with Us = , we have for a smooth gauge field 


5W[C:] = -z96a{A:W[C:] - VF[C:]AD. 


(34) 


This gives (if the eigenvalues of the Wilson Loop are not degenerate) 


tan 2/3C 


e2*7r 


2dag^-4{A^sW[C^s] - (A;W[Q] - W[C-sWs)i;2 

"Ai - A 2 - igSa4iA-,W[Cl] - IF[C r]A;)V;i + 4(A^IF[C;] - W[Cl]Al)i:2igSa 

JlSag^J (Ai - A2)Vj(Aj)V’iV’I^sV'2 

Ai — A 2 

=sign(Ai - X2)5ag^JtT:[{T^ + /T2)6lt Aj6»]tr[(ri - AlO], 

tr[(Ti -h A^6l] 

“tr[(ri 


/3 is proportional to 5a, so, baring changes in the phase, the eigenvectors evolve smoothly around 
the Wilson Line. 
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Considering changes perpendicular to the Wilson Line, we have 


= tgSaSr V u:u:^s. ■ • • driu:^s.)iu:+^sJ ■ • ■ iu:^sJm^w[c:] + w^], (se) 


which gives, since A 1 A 2 = 1 for W S SU{2), 


tan 2/3’’ =2 


5cj5rg^Y.s',s" e"'*tr[(ri + iT^)e\,drAs.es']tr[{T^ - iT^)6\„drAs„es"] 


Ai - A 2 + igSaSrj:^, 
SaSre-^^f:' tT[{T^ + tT^)el,drA,,es'] 

T,s' SaSre^^^"' '^'*"'“^"tr[(ri - iT'^)el,drAs>es'] 


(37) 


We see that, expect when the eigenvalues are near degenerate, /3 will be proportional to Sr — 
though not 6a as s" 0(1)- This means that the 9 field is discontinuous when the 

Wilson Loop has degenerate eigenvalues - and we shall investigate the consequences of this later. 
For now, we just consider the other situations with infinitesimal /3, and ask if it is possible that 
there can be a large change in 9 without a large change in (3. We obtain 


9r+Sr 

Tj _ f COS/S’" — sin/3’’e“*'>'’^’*’ 

sin/3’'e*'’'’^’*’ cos/3’' 




(38) 


where the phase S is chosen so that whatever fixing condition we have chosen is satisfied. We 
denote as the appropriate transformation when performing a small translation in the direction 
of increasing s and when translating in the direction of increasing r. In section [5l we shall 
use these relations in our discussion of string breaking. 

There is one further discontinuity as we evolve 9 over s and r. We may parametrise 9 as 


/ cos a fsinae’’’ A / ^ 

y isinae”*'’ cos a / V ^ ) ’ 


(39) 


with 0 < a < 7 r/ 2 , c G R and d G R. Note that the parameters a, c and d are not gauge invariant, 
so this parametrisation is only defined after we have fixed to a suitable gauge. The phase d makes 
no physical contribution, and should be fixed by some suitable condition, but for completeness we 
still give it explicitly. We perform a small translation in either the r or s direction, and obtain 9' 
which is parametrised by a', d and d! 

cos a cos /3e®‘^ + i sin a sin 
cos a sin + i sin a cos 



— cos a sin/3e*^‘^ + i sin a cos/3e*^’’ A 

cosacos/3e“’'^ — f sinasin/3e“’*^’’+'''“‘’*^ ’ 

(40) 


which allows us to read off the new parameters o', d and d\ 


cos a! 
tan((i' — d — S) 

tan(c' + d — d' + S) 



2 cos a cos fi sin a sin /3 sin(c + 7 — 2d) 


sin a sin /3 cos(c + 7 — 2d) 
cos a cos /3 — sin o sin /3 sin(c + 7 — 2d) 
cos a sin /3 cos(2d — 7 ) + sin a cos /3 sin(c) 
— cos a sin (3 sin(2d — 7 ) + sin a cos /3 cos(c) 


(41) 


It is clear that when cos a and sin a /3, we have a' = a + 0{(3), d = c + 0(/3) and d' = d + 0(/3), 
which means that 9 is differentiable with respect to the spatial coordinates at least once. However, 
if a ^ 7 r /2 we have d' = c + j — d + S+{i'+ ^) 7 r and c' = 2c — 2d + 7 + (i^ + \)tt, where v is an 
integer, and, at a ^ 0, c' ^ 2 d — 7 + (z/ + i) 7 r. At both a ^ 0 (where 9 is diagonal) and a ^ 7 r /2 
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(where 6 is off-diagonal), the parameters which describe 9 are discontinuous even when /3 is small. 
a itself, however, remains continuous and the eigenvalues of 6 are dependent only on a once we 
fix the phase d. 5 can be chosen according to our rule for fixing the U(l) phase d. For example, 
if we enforce d = d' = 0, then 6 = — arctan ( - sinasin/ 3 cos(c+ 7 ) \ ^ 

’ y COS a COS p— Sin asm p sin.(c+ 7 ) J 

This discontinuity of 6 may be in principle be at a single point, or there may lines in the 
surface spanned by r and s where these discontinuities occur (if they occur), and these lines 
will either be open, form closed loops or will be infinite in extent. If they are lines, then we 
require that a' = a = 7 r /2 (for example), which means that 9 will be unchanged as we make a 
small displacement in space in some arbitrary direction. If 9 is unchanged, then the Wilson Loop 
must also be unchanged as we move across space time. However, baring a remarkable gauge field 
configuration, it is unlikely that this will occur for any more than an infinitesimal distance. We 
may as well then neglect the possibility that these discontinuities are in lines, and focus on the 
case that they are at isolated points. 


3.2. SU(3) 

In SU(3), we introduce the following operators which form a complete basis of the Hermitian 
traceless generators of SU(3) 




<('2 






0 

0 

0 0 

0 0 

0 0 

0 

0 0 

0 0 

0 

1 0 0 

0-10 

0 0 0 



(42) 


Any SU(3) matrix, 9, can then be parametrised in terms of the eight quantities oi, 02 , 03 , ci, C 2 , 
C 3 , d^ and ds as 


Q _gja303 gio202 giai <pi 

As in SU(2), these parameters are not gauge invariant. 
Multiplying out the various components of 9 gives 


/ 011 

012 

0l3 \ 


0 

0 ^ 

021 

022 

^23 

0 

id^—ids 

e n /3 

0 

V ^31 

032 

^33 / 

1 0 

0 

e / 


with 

011 = cos 02 COSOi 

012 =zcoso2sinaie*°^ 

013 =z sin 026 *°^ 

021 = —cosoi sin 03 sin 026 **'°'*“°^^ + z cosos sinoie 
022 =coso 3 COSOi — z sinoi sina 2 sino 3 e*^'^^'''°'®“°^^ 
023 =Z sin 03 cos 026 *'^'* 


(43) 


(44) 
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031 =* cos as cos ai sin 026 — sin 03 sin me *('=i+‘=3) 

032 =zsin03 cos016“*'^^ —cosa3 sinoi sin026*^'^^“'^^^ 

^33 = cos 03 cos 02 ( 45 ) 

Note that ds and ds can easily be constructed from the phase of the 11 and 33 components of 
0 , 62 and 02 from the 13 component once the phase has been removed, and 03 and 63 and oi and 
d from the 23 and 12 components of 0 . We can again construct the differential of 0 with respect 
to a coordinate by multiplying 0 by a matrix B constructed in the same parametrisation as 0 : 


/ / 3 ii 

l 3 12 

Pl 3 \ 

/ 


0 

0 ^ 

B=\ 1321 

1322 

1323 


0 


0 

\ (331 

1332 

^33 / 

V 

0 

0 

6 Vs y 


with Pij taken from equation ( 1451 ) but with the substitution 6 ij —>■ jdij, ai ^ ( 3 i, Ci ^ ji. 

Again, except in the neighbourhood of degenerate eigenvalues of the Wilson Loop, /3 will be 
infinitesimal and we can neglect terms of 0 (/ 3 ^) and higher. With the shifted 0 ' = OB, we have 


011 = (cosoi cosa26 v /3 ^ cosa2 smai6 ^ —,82 sin 026 ^ 1 

rf ( ■ Q d 7 l + '^+^ 3 ) ■ • dci + -^—^3) Q . hc2—73— 

012 = UCOSOI COSa2/3l6 ' + ^C0S02Smal6 —/33Sm026 ^ 

0(3 = ^icosai —/33 cos 02 + j sina2e*^‘^^~^^ 

023 = cos 03 sinoic”*'^^ — cosoi sin 03 sin 026*^'^^“'^^^)+ 

1/336*^"*'^''’^”^^^^(cos03 cosoi — zsinoi sin02 sina3e*^'^^'''°^“'^^^)+ 
i sm 03 cos 026 ^ W/a M 6 


2(53+Z—^ 

e ^ 


eVs 


e 


^33 COSOs cosoi sin026 — sinossinoie dc3+ci)^_|_ 

iP3e^^~^^~^^~‘^^\isma3 cos 016“*'^^ — sinoi sino2 cos 036*^'^^“'^^^)+ 

_2?48_ , , 

cos 03 cos 026 16 v^, ( 47 ) 


where, for the d = 0 fixing condition, 6 is tuned to ensure that 0 (i and 033 are real. Once again, 
we see that o^, Ci and di will be differentiable at least once where /3 is continuous except where 
either oi = 7 r/ 2 , 02 = 7r/2, 03 = 7r/2, oi = 0 , 02 = 0 , or 03 = 0 , where the parameters Ci are 
discontinuous (possibly because the 5 required to fix d is not infinitesimal). 


4. Topological solutions 

4 . 1 . SU( 2 ) 

Now suppose that contains a non-analyticity. We integrate the field around a loop C 
parametrised by a surrounding the discontinuity, bounding the surface integral by an additional 
line integral f^ddu^-, far enough away from the discontinuity that u^~ is analytic along C. We 
define {Cn} as the set of curves surrounding all these discontinuities, and S the area bound within 
these curves. We can write 




x-\-5a 5 


( 48 ) 
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and since u is continuous on C, after fixing the gauge we can expand U = 1 — i^gSd-A°‘X°‘ and 
9l.9j;+sa = l+5a9^da9. We define Xq = ^9^{X+X^)9. For example, in SU(2) we may parametrise 
X as 


1 = 


cos a; zsinxe**^ 
fsinxe”®*' cos a; 



(49) 


then X + It = 2 cos a: cos ui/ = 2Xo, where I is the identity operator (this expression will obvi¬ 
ously not be exact in higher gauge groups). If everything is analytic along Cg, and we consider 
infinitesimal displacements, we expect X and U to be close to the identity operator, and x and w 
will both be O(dcr). Xq is thus I -I- 0{5a'^). 

This gives 


i5dul, = ( tr [X9lXljJlU^,,9,+saf. 


tr(AJ)2 

^^-^tr[A^dt (It _ x)9, _ i^A^dddt (Xt + X)gX;_,A“d, + A^0t (^t + X)9Ja9ld,9]. (50) 


Using (HU the first term in equation (I50p gives zero, while the second term will not contribute 
to an integral around C \iU and Xq are continuous and the area of the loop is small enough. We 
therefore concentrate on the contribution from the final term. Equation (1201) is then replaced by 


9lW[Cs\9s 




(xes)n(x^i:) 




(51) 


where dT, is an element of area. 

In SU(2), after fixing to a suitable gauge (where the gauge field U is smooth), we parametrise 
9 as in equation dMD, 


9 =(cos a -I- i sin 



0 


(j) = 





(52) 


with cGK, 0 <a<7 r/2 and determined by the fixing condition. As this contains the 
eigenvectors of IU[C's], it is differentiable except where IU[C's] has degenerate eigenvalues and 
those locations where a = 0 or a ~ 7r/2, where c is ill-defined. The parameter c may wind itself 
around those points where a = 7r/2ora = 0. We may parametrise the plane of the Wilson Loop 
in polar coordinates (r, </)), with the origin at the location where a = 7r/2ora = 0. The continuity 
of 9 away from r ^ 0 implies that at infinitesimal but non-zero r, c(r, </) = 0) = c(r, (jj = 27r) -|-27rF„ 
for integer winding number and if c is ill defined at r = 0 then we may find that Vn ^ 0. 

It is straight-forward to calculate 

9’^drj9 = [idcracj) + iX^d^ds -I- i sin a cos a4>daC — i sin^ aSacA^] . (53) 

We integrate around a curve at fixed a, and we may choose a fixing condition which keeps d 
constant. This leads to 

9^W[Cs\9s =6*"^^ k-n do-sin'^ ad^c^ti[Xon]^ 

[/(„gs)n(xg£) 27rsin2 ai/„A3] ^ 

where the last equality is only valid if Xq is close to the identity operator, as we expect. In SU(2) 
(but not higher gauge groups, where the expression depends on more parameters) we discover 
that were we to use an infinitesimal loop around the discontinuity, then the area integral gives 
9\W[Cs\9s ^ g27rzj/A Jqj. integer v, which cannot contribute. However, if we integrate at a larger 


13 







radius, where a has fallen to some average background value oq, we do see a contribution to the 
Wilson Loop from the discontinuity, proportional to 27rr'„ sin^ oq, and we may still find a large 
contribution to the string tension. The structure in the field is extended over a small region of 
space rather than being a (5-function spike. This conclusion may also be reached by considering 
gauge invariance: while -F is gauge invariant, 0 is not and therefore the precise location of the 
discontinuity in 6 will depend on the gauge (although it should not depend too strongly on the 
gauge). The peak in F cannot therefore be precisely at the discontinuity in 9 in every gauge, so it 
cannot be a i5-function peak around that centre. We shall discuss gauge invariance in more detail 
in section HU 

We will refer these objects in the field strength as CDG merinthons. They are topological 
objects in the homotopy group tti (distinct from monopoles which have the homotopy group 7r2, 
and instantons which are of the group tts). They can be identified by a large value of T", a = 7r/2 
or a = 0, and a non-zero winding of the parameter c around the merinthon. 

It is also important to note that the merinthons at a = 0 and a = 7r/2 with the same winding 
contribute with opposite signs. If we compare the results for the integrals of the restricted potential 
around fixed curves at a = (5 and a = 7r/2 — 5, for some small 5, we find that the term proportional 
to the winding number gives 



gAfj,dXfj, =A3 sin^ 62Triy -!-... = \^5‘^2'kv -|- ... 

gA^dx^ =A 3 sin^(7r/2 — 5)2'itv -!-... = A3(l — 5‘^)2ttv + ... 


( 55 ) 

When this is multiplied by i and exponentiated, = 1 can’t contribute to the Wilson Loop, 
so we can neglect this part of the expression and write, 


® gAfj^dXfj, = —Xs6^2tti' -|- ... ( 56 ) 

J a—7Tl2—6 

Thus with the same winding number, the two topological objects contribute with opposite signs 
to the Wilson Loop. This is important as we consider how the merinthons with positive and 
negative winding numbers can be distributed across the lattice (see figure [IJ . Assuming that 
c is continuous apart from at the merinthons, then we must have a positive winding number 
object next to a negative winding number object. Thus two a = 7r/2 objects next to each other 
(with opposite winding numbers) will give contributions to the Wilson Loop which approximately 
cancel (ultimately leading to a total contribution proportional to the square root of the number 
of objects, i.e. a perimeter law scaling for the Wilson Loop), and the same can be said for two 
0 = 0 objects neighbouring each other. However, if we have a, a = n/2 merinthon with positive 
winding number situated as a spatial neighbour to an a = 0 merinthon with negative winding 
number, then the contributions of the two objects would be additive (ultimately leading to a total 
contribution proportional to the number of objects, i.e. an area law scaling)!^ 

The total Wilson Loop will be the product of a boundary term and contributions from all the 
merinthons contained within the planar Wilson Loop. As we can expect the number of merinthons 
to be proportional to the area of the loop, this leads to an area law scaling for the Wilson Loop 
and a linear string tension. 


4.2. SU3 

The SU(3) system is considerably more cumbersome than the SU(2) theory. We start by 
constructing the operators = 9X^9\ n® = 0A®0^ and O^d^O, where we parametrise 9 = 


®An alternative, which we don’t consider in this discussion, as that there is another type of topological object, 
such as a monopole, neighbouring the merinthon. 
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Figure 1: An Illustration of how topological objects can be distributed across a plane. The cross indicates the 
centre of the topological object. The circles and arrows indicate the winding of the parameter c around each 
topological object. In this simplified picture, it is clear that if c is continuous outside the topological object, the 
arrows of neighbouring objects must point in the same direction at the point of contact between them. In this 
simple picture, each object of positive winding number (clockwise arrows) must have objects with negative winding 
number neighbouring to it. While this can be avoided by having the direction of flow of c doubling back on itself, 
such a scenario would require that the flow should be three times as large on another part of the circle, meaning 
that the problem is not eliminated but just transferred elsewhere. 


g*a303e“202gmi0igid3A3+zd8A«^ ’yyg 


9^d^e =Si^ + ^2^ + S 3 ^ + 54^ 

C id 3 \^+ids\^ 

_g-id3A^-id8A*g-iai<^i^^^giai<^i^gid3A^+id8A® 

—^-id3X^-idgX^ ^-ia24>2 Q ^gia2</'2 ^giai<^'l gids A^+idsA® 

_g-id3A^-id8A*g-iai<^ig-ia2 02g-ia3'^3^^^gia3'^3^g*a2 02g*ai<^'lgid3A^+id8A® 


where explicit expressions for vr and , ^ 4 ^ are given in Appendix C 

Close to the merinthon we can write the Abelian decomposition as 

TT _ V fl iSc!X^U3^+i&(7\^U3a 

Ua — ^crOae 

which gives 

gZ5<TA3u3„+i5aA»fl8. = 9 lX^U0a+Sa. 


(57) 

(58) 

(59) 

(60) 
(61) 


(62) 


(63) 


The two components of the field strength will be given by the integral of M 3 or ug, around a loop 
surrounding the merinthon at an Oi that have fallen to some background level. If the fields X and 
U are smooth, and with 9 differentiable along this path, we have 


uga = - 


(64) 


which implies that (using equations (IC.4|) - (IC.8|) 1 


ul = dtjdg — sin^ Oi^^ci + 


sm 02 


COs(2oi)a<^C2 + dcyCg 


sm 03 


( 003 ( 202 ) — 2)cos(2ai) ) — 


drjCg ^ ^ sin 02 sin 2 oi sin(ci + C 3 — C 2 ) ) + d^jag sina 2 sin( 2 oi) cos(ci + C 3 — C 2 ) 


Ufj — I dfjdg d^C 2 


•\/3sin^ 02 


2 + cos(2o2) 1 . 


(65) 
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The components of the field strength tensor are then given by (excluding terms which will not 
contribute either because they are multiples of 2 tt or zero) 

r / -2 -2 

= j>_da(- sin^ aid^ci + cos(2ai)9^C2 + (cos( 2 a 2 ) - 2) cos(2ai)- 

sin 2 a 3 , • n • / , \ 

o^j.C 2 — 2 — sma 2 sm 2 ai sm(ci + C 3 — C 2 ) 

= j) da ^-i9^C2 cos( 2 a 2 )^ . (66) 

where dS is an element of area which covers the merinthon and C the curve which surrounds the 
merinthon. Only those c^s which wind around the merinthon will vary on this infinitesimal loop; 
all the other parameters can be treated as constant. These field strengths are non-zero if one 
of the Ci has a non-zero winding number around the merinthon, and they will not always be an 
integer multiple of 27r. In particular, for the terms proportional to d^C 2 and we may pick 
infinitesimal loops where the a parameters do not vary significantly. The merinthons thus will 
lead to a non trivial structure in the restricted field strength, and if the number of merinthons 
within an area is proportional to that area, an area law for the expectation value of the Wilson 
loop and confinement. 


4-3. The Electromagnetic field 

We saw in equation (l20l) that the field strength can be written as, 

- d^ul, (67) 

with (from equation (lA.171) 1 

uj, = -^tT{n^igA^-i0dfj,0^)). ( 68 ) 

In 0, the potential is split into two parts Ae = tr A^ and Am = 'n^0d^0fi The field 

strength, 

Pfiu = ^ [9^tr(n'’A^) - dviT:^A^)] + (69) 

is often split into two components]^ 

=\ [9^tr(n'’Ay) - dAT:WA^)\ 

='^tr(n'^[5^n'‘,a^n'']). (70) 

The magnetic current can be defined through Maxwell’s equation 

(71) 

and while d^Ep^ = 0, as in electromagnetism, in general 8,^11 pa ^ 0. This means 

that this field strength, coming from the colour field rather than the gauge field, is associated with 


^In the literature, E is known as ‘electric’ and H as ‘magnetic’, but we will not be following this terminology 
but refer to the ti components of the fields as electric and the ij components as magnetic, where i and j are spatial 
directions. For a monopole, H is purely magnetic and the terminology is sensible; however merinthons are not 
monopoles and will usually also carry an electric field. 
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a magnetic source. Note, however, that it is not in general a pure monopole (except for certain 
specific choices of 0), since the electric current, defined through d^Hf^u ^ 0. Neither nor 
are by themselves gauge invariant, but only F = E + H. We expect that close to the merinthon 
solutions F will be dominated by H since this is the part of F which is sensitive to the winding 
of c around the centre of the merinthon. In analogy to QED, we shall refer to the components of 
F (and also H) as the electromagnetic fields. 


4 . 3 . 1 . SU(2) 

In SU(2), we see that (using equations (IMl) and (IA.25I) and setting d = 0), 


i 

^9 

sin(2o) 


= — — tr [dfj,a(j) + sin a cos — sin^ , 9^a(/) + sin a cos acjid^c — sin^ a9^cA^] 


9 


{d^ad„c — d^cdyo). 


(72) 


The Yang-Wu monopole 2^ 26, 271 is one classic magnetic monopole solution for winding number 
z/ = 1. Here we parametrise the coordinates as (with the origin at the location of the monopole) 


( t \ 


/ cos 03 \ 

X 


sin 03 cos 02 

y 

^ T 

sin 03 sin 02 cos 0i 

V ^ 


y sin 03 sin 02 sin 0i J 


(73) 


The reason we choose this parametrisation is so that we can easily respect the symmetries of the 
theory. Since we are considering point-like objects, it is most natural to use a spherical polar 
coordinate system. QCD, after averaging over the gauge links, possesses a spherical symmetry, 
which we have broken by inserting the Wilson Loop in the xt plane. However, we should still be 
free to rotate the coordinate axes in the yz plane. This particular parametrisation allows us to 
absorb such rotations into a shift of ^1. If, however, we parametrised the coordinates differently 
(for example by exchanging t and z in equation (17^ 1. a rotation in the yz plane would require 
some non-trivial change of all the angular variables. The only changes we can easily make to 
the coordinate parametrisation are to interchange the t and x axes or rotate in the yz plane (we 
shall only consider rotations of 7r/2 since we are primarily interested in comparing with a lattice 
theory). 

For the Yang-Wu monopole we select the solution with c = tpi and a = '02 /20 This leads to the 
solution for the electric and magnetic fields E = 0, B = where rg = {x, y, z); in other words 

the true magnetic analogue of a point like electric charge. Note that there is no time dependence 
in the monopole field strength: these objects will reveal themselves as stings of electromagnetic 
flux extending in the time direction. 

The literature concentrates on the possibility of finding these monopoles in the colour field, 
which are then assumed to condense and provide confinement through a dual Meissner effect. 
We will later, therefore, in addition to the merinthon solution suggested above, consider what 
structures may appear in the fields if these Yang-Wu monopoles dominate. 

In our case, we have something a little different. We will start by considering the coordinate 
system given in equation (1731) . We consider the objects close to a = 7r/2 here; the objects at a = 0 
can be treated in a similar way. We can then write the variables, at all locations away from the 
singularity, as a = 7r/2 — G'(r,-01,02,'03) where G is zero at r = 0 and c = i?(r,0i,02,03); and 
for a winding number one merinthon we require i?(03 -I- 27r) = 77(03) + 27r. The functions R and 
G will depend on the gauge, but in a way which allows F^ to remain gauge-invariant. This gives 


®Here one would have to choose the coordinates parametrising 9 within the ranges 0 < a < tt, —- 71/2 < d < 77 / 2 , 
c S R, rather than as we set it out earlier. Note that for this monopole, 9 is discontinuous (0('!/’2 = Stt) = —9{'Lj)2 = 0)) 
but V? remains continuous. 
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US 


tan'll =- 


y 


tan'(/;2 = 




tani/ja = 


^/x‘^ +y^ + z"^ 


=\/ a;^ + 2/2 + 


(74) 


and 


9V'i cos 'i/'i y 

dz r sin '02 sin 03 
901 sin 01 z 


9?/ r sin 02 sin 03 z^ + ?/^ 

902 sin 01 cos 02 a;z 

9z r sin 03 y^z^ + j/2(a;2 ^2 

902 COS 01 cos 02 xy 

dy r sin 03 ^z^ + y2(a;2 + y2 ^2) 

902 sin 02 \/ + y"^ 

dx r sin 03 a;2 + i/2 _|_ ^2 

903 cos 03 sin 02 sin 0i 
dz 


zt 


y^a;2 _|_ y2 _|_ 22(a;2 _|_ y2 _j_ 2;2 _|_ ^2^ 

903 cos 03 sin 02 cos 0i ?/t 

9?/ r 

903 _ cos 03 cos 02 
9a; r 


^ a ;2 + ?/2 _|_ 2 ; 2 ( a ;2 + 1/2 _|_ ^2 -)- 0 ^ 

xt 

Y^a;2 + y2 _|_ 2:2(a:2 + j/2 + z2 + 


903 sin tfj^ ^a;2 + 1/2 _|_ 2;2 

9i r (a:2 + i/2 _|_ ^2 + /2)' 

The simplest possibility is to set G = G{kr) for some k and i? = i?(03), neglecting the 0i and 02 
dependence. We shall consider a more general solution later. Using equation (ED), we obtain, 


1 


= Ei =- sin(2G)i?'G' 


,, , / 9r 903 dr dips 


dz dt dt dz 


sinf2G) 

:- -R'G' (— sin^ 03 sin02 sin0i — cos^ 03 sin02 sin0i) 

sin(2G) ^,^, z 


(75) 


\/a;2 + 1/2 _|_ ^2 

where i?' and G' indicate the first derivative of these functions and is the electric field. Similarly, 

y 


3_£^3_ 


TJO _ TT'O _ 


9 r 


^3 ^3 1 sin(2G) 

—^x — ^ ^ 


9 r 


V ^ a ;2 + 1/2 + 2;2 
X 

\Jx^ + y2 + z2 


(76) 


For the magnetic field we find. 


= -Bt =-sin2GG'i?' 
^ 5 


1 __ I of f z dips _ y dips 

r dy r dz 


1 sin2G , , 


9 r 


G'R' (sin 03 cos 03 sin^ 02 sin ipi cos ipi — sin ips cos ips sin^ 02 sin0i cos 0i) 
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=0 


Hi, = -Bl =0 

Hi, = -Bl =0. (77) 


Here, the merinthon generates a pure electric field. This is, of course, not the only way we can 
parametrise the coordinates. For example, we could have chosen the coordinate system 


(G 


/ sin ip 3 sin ^/>2 sin ipi \ 

y 


sin ip 3 sin ip 2 cos ipi 

t 

^ r 

sin 'ip 3 cos ip 2 

\x ) 


\ cost's / 


which in effect swaps t and x. This gives a different solution for the field strength tensor 


fjd _ fp 6 

^xt 


Tj^ d3 

tiyx - -^z 

ttS d3 

^zx - 


1 sin(2G) t 

9 r ^y2_^_^2 

_ l sin(2G) ^,^, y 

9 ^^2 _|_ y2 _|_ ^2 

_ 1 sin(2G’) z 

9 r y2 y ^2 ’ 


(79) 


with the other components of H^ zero. We could also use any other rotation of the coordinates 
consistent with the space time symmetries. 

For a general function G{r,ipi,ip 2 ,'ip 3 ) and R{'ipi,'ip2,'4’3) (neglecting any radial dependence of 
c, which may in practice only be valid for small enough distances and certain choices of gauge), 
with di = djdipi and = y'^ + and + z^, we find for the parametrisation (1751) 


Hl^ = - 1 sin2G( drGd^R— - {diGd^R - d^GdiR) 
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(80) 


We must also consider every permutation and rotation of the coordinates to obtain other possible 
solutions. 
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At small r, each of these terms gives a finite and non-zero result if the derivatives of G and 
R are non-zero and both G and d^^G are proportional to r. Our ultimate goal is to measure the 
electromagnetic fields on the lattice to try to identify these structures, and this means evaluating 
how the fields behave at large and intermediate distances. This is complicated by the unknown 
functions G and R (which are also not gauge invariant and therefore uniquely defined). We require 
that at least one of the dip^R ^ 0 (and, in practice, on average take an integer value). We can, 
however, study the explicit coordinate dependence within equation (ISHll for G ~ r (obviously this is 
in practice only an approximation, as we only consider part of the expression for the field strength, 
but nonetheless we may expect these features to be visible to a certain extent). We can see three 
types of behaviour at large distance. Firstly, there are terms such as {diGd 2 R — d 2 RdiG )-—^— 

^yz^xyz 

(taken from Hy^) which fall away at least as 1 /zi for large spatial distances, while increasing 
at large t if diG is linear in r - we find a possible solution where the field grows with distance, 
although we may hope that these are excluded in practice. Secondly, we have terms which form a 
line of large electromagnetic field in one of the directions, for example in we find the term 


-sm2GdrGdiR^ <-drGdiR^, (81) 

9 rr^z 9 rr^z 


and the maxima of this function fall off at large distances from the origin as 1/z^, 1 /y^ and 1 /t 
but remain roughly independent of x. This will correspond to a line of large magnetic field (which 
may be broken or change sign) in the x-direction, and we will label it as a /i-string, where y 
indicates the direction (in practice, the line of field strength will only have a finite extent because 
of the G and R dependence). We do not mean anything by this terminology beyond that the 
large field may extend in the /r-direction and falls of rapidly in the other directions. We require 
that the leading order r dependence in G should be linear if we wish to avoid zero or infinite field 
strengths at r = 0. Alternatively, by permuting the coordinates, we also find these strings may 
exist in the other electromagnetic fields. Finally, we have what we will label as points where the 
electromagnetic flux falls as at least 1/r in each direction at large distances, r > 1/fc, for example 
the term sm 2 GdrGdiR^ in Hy^. 

We will concentrate on solutions with diR = 0 and 82 R = 0. Our reason for doing so is that 
we expect the discontinuity to occur in a point; thus there is no obvious reason why the objects 
should wind around this discontinuity in the yz plane and not the xt plane; if the origin is at the 
discontinuity and c proportional to the angle we might expect c[— 6 xy) = c{5xy) + tt, and only 
setting c = '03 accomplishes this. If there is no winding number related to variations in ijji and ■02, 
then we might expect any fluctuations in the fields in these directions to cancel out as we integrate 
around a loop surrounding the discontinuity. Also, only winding in the xt plane can contribute 
to confinement. In our numerical simulations, we only use Wilson Loops in the xt-plane, so we 
restrict ourselves to these loops here. This means that we should expect a symmetry between the 
structures in the x and t components of the electro-magnetic field, and also between the y and z 
components of the field, but not necessarily between (for example) the x and y components. This 
all suggests that it is safe to neglect this angular dependence in R. With this simplification, we 
obtain. 
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Equation (|73)) 

point 

x-string 

x-string 

t-string 

t-string 

t-string 

t X 

point 

x-string 

x-string 

x-string 

t-string 

t-string 

y^z 

point 

x-string 

x-string 

t-string 

t-string 

t-string 

t O x,y -O' z 

point 

x-string 

x-string 

x-string 

t-string 

t-string 


Table 1: Some Possible configurations of the restricted electromagnetic fields; every time we specify a string there 
may also additionally be points. This table lists the fields which could in principle be present; it is not necessary 
that each of them would be visible for a given merinthonin practice. 
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t-string 


1 


= Bl =-sin2G( diGd^R 
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H^^=BI = -- sin2G( diGd^R 


zxt 


/j^2 .y’2 v* 

' 'yz'xyz 


t-string 


(82) 


In this parametrisation of the coordinates, the electromagnetic fields will therefore appear either 
as points (not necessarily spherically symmetric, but falling in every direction) or strings in either 
the X or t direction. We also consider the three other parametrisations of the coordinates which 
arise from the exchanges x ^ t and y ^ z. We summarise the results in table [H With these 
choices of coordinate parametrisations, we see that Ex will only have points, Ey and E^ strings 
along the x axis, Bx may have a string along the x or t axis, while By and B^ strings along the 
time axis. Of course, in practice this picture will be obscured by fluctuations in the fields and 
effects from the unknown functions G and i?; we also cannot say how long the strings will be. 
It also depends on that G has some angular dependence, and in particular has some variation in 
the yz plane (as parametrised by the angle V'l)- Equally, the strings depend on different angular 
differentials of G; if G has no angular dependence they will not be present, and if it depends on 
'tjj 2 but not on tpi only some of these strings will appear. We will search for these patterns of 
electromagnetic fields in in lattice QCD in section 

We can also consider the patterns we can expect from the Wang-Yu monopole and other 
solutions for the fields. We have already seen that the monopole solution with the coordinate 
parametrisation given in equation (1731) gives t-strings in the magnetic field; we may also study 
the other options when we permute the coordinates. Additionally, we list a few similar simple 
alternative forms for the functions R and G. These are summarised in table [2l It should be 
observed that if a is linked to an angular variable then either Ex = 0 or Bx = 0; while if c is 
linked to the angular variable '4’2 rather than t/js with G ^ r, then the Electric fields have t-strings 
and the magnetic fields x-strings: the opposite of our expected model in tabled] We note that a 
combination of a = ' 03 / 2 ; c = V'l monopoles and a = /2] c = 'ip 2 structures may have the same 
pattern of points and strings as the conhguration of table dl the difference is that there is in this 
case no correlation between the spatial locations of the strings and the points in the E^ field. 

That the structures in the field strength are not point like but contain vortices is important 
because it means that we cannot evade the large field strength by distorting the surface E. For 
example, we may create a bump in E to avoid the merinthon and the structure in the Ext held: 
however, this will require integrating over a surface in (for example) the yt plane where Eyt is 
large, and the yx plane where Fyx is large. This means that it is not naively ruled out that the 
surface integral bound by the Wilson loop is not affected by the choice of surface to integrate over. 
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Equation ( 1751 ) 
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Equation ( 175 )) 
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t-string 
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point 
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x-string 
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x-string 
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t-string 
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t 0 x,y 0 z 
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point 
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Table 2: Some Possible configurations of the restricted electromagnetic fields for the parametrisations a = t^2/2, 
c = 'ipi (including the Wang-Yu monopole) (top); a = -03/2; c = -02 (second top); a = -03/2; c = 0i (second bottom); 
G = kr; c = -02 (bottom) 


4.3.2. SU(3) 

SU(3) gives a similar picture as SU(2); here we only outline the argument about which struc¬ 
tures will emerge in the field strength. The starting point is equation (16411 . and by differentiating 
this we can construct the two field strengths and F®. For example, we find that 

= 9),ul - d^ul = -V3sin{2a2)idfj,a2d^C2 - d^C2d^a2)+ 

-y3cos(2a2) sin(2a3)(9^a3i9jyC3 - d^C3d„az)+ 

A/Ssin^ 03 sin(2a2)(i9i/C3i9^a2 - d^C3dua2). (83) 

The first two terms are similar to the quantities already discussed in SU(2), and the same analysis 
applies. The third term, mixing the differential of C 3 with that of 02 , is new and we need to 
consider this as well. is similar to F®^,, with similar expressions arising, plus additionally terms 
proportional to the product of the differential of to of the qs, for example 9 ^ci 9 ^C 2 . We therefore 
expect to see the same structures as in SU(2), plus possibly some additional structures from these 
extra terms. For example, in 0® with C 3 ^ '03, we find the contribution ^^10 ^ 

while in F®i, we find the terms - 2 ^* . In general 02 , unlike 03 , will not be at 

its minimum or maximum value at the merinthon . There is therefore no reason to require that 
002 / 501.2 will be proportional to r at small r (we can discount the possibility that C 2 and C 3 both 
have winding around the same point). We can also neglect terms such as d^Cid^C 2 because we do 
not expect both ci and C 2 to wind around the same point. This means that these additional objects 
in the electromagnetic field strength will just be points. The same argument applies to all the 
additional terms in the F®^,. Therefore, the electromagnetic structure in SU(3) will be the same 
as in SU(2) but there might be additional point-like structures appearing in all the components 
of the electromagnetic field; this applies equally for both of the electromagnetic field tensors, F®^ 
and F®^. 

4 . 4 . Gauge Invariance 

is gauge invariant, which is what we expect for the mechanism behind confinement. However 
9, X and U are dependent on the gauge. This presents a conceptual difficulty, since we have stated 
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that the gauge invariant field strength is influenced by objects which arise due to features in the 
gauge dependent parametrisation of 6 . The key contribution, however, comes from the winding of 
c around the merinthon, and we here show that this is gauge-invariant (at least in SU(2)). 

Since this formalism requires that the gauge field is continuous, we can only use continuous 
gauge transformations. However, to undo the winding in 9 we require a discontinuous gauge trans¬ 
formation. Thus the discontinuity in d at a = 7 r /2 will survive any smooth gauge transformation, 
and the merinthons will remain. For example, in SU(2), we can parametrise an infinitesimal gauge 
transformation as 


A = 


cosli isin^ie*^^ 
tsinZie”**^ cosh 


0 \ 

0 J ’ 


(84) 


with h and h infinitesimal and 0 < ?2 < 27r. Recall that under a gauge transformation 6 —>■ Ad. 
Multiplying Ad together (choosing = 0) gives 

. „ _ / cos/i cosa — Zi isinae*^'^''"^*^)-I-zcosaZie®*^ \ / e**'® 0 \ 

z sinae“*(‘^+^'=®^-I-z cosa^ie”*'^ cos/i cosa —/i sinae“*^'^“°“^'='^ / \ 0 )' 

(85) 


If we define 


cos a 



2li sinacosacos(l2 


c) -I- /f sin^ a, 


then we can rewrite this as 


( 86 ) 


Ad = 


cos a 


zsma e 
cos a' 


COS li cos a—li sin ae 


i{l 2 -C-2Z3) 


COS/i COS g—/i sin ae 
cos a' 


Sin a e 


(sinae*(‘^+^'3^-I-cosa^ie*'=)(cosa — sinae*('= ° 
cos a' 


(87) 


We fix ds = 0, so we remove the matrix on the right by a 17(1) rotation. From the remaining 
term, we can read off as the phase of the top right component of the matrix on the left. In 
particular, we see that 


/ i(c'-c) _( 


sma'e"'" =(sinae*^^*^^-I-cosalie*^*^ °^) 


cosa — b sinae*^^^ ° 


cos a 


, , , smacosa , , cos^ a ,, , , sin^ asin(l2 — c) , , 

sin c' - c) -- sin213 + h- -;-7 sin(l2 - c) - li-^- / , . 88) 


sm a' cos a' 


sm a' cos a' 


sm a' cos a' 


If a -K 12 and a 7 ^ 0 , we can read off, to lowest order in li and I 3 and given that c = c' at 

li = I 3 = 0 (i.e. there is no gauge transformation), 


a' — a =li cos(l 2 — c) 

c' — c = 2/3 -I- h sin(l 2 — c) cot a — /i sin (/2 — c) tan a (89) 


The winding around the merinthon becomes A^ 5 tr[Aro] j> d^c'dd = A^ 5 tr[ 7 fo] j> 9a (c' — c)dd' -I- 
A^TTz/trXo, and since li and I 3 are continuous functions and c only contributes to (c' — c) within 
a trigonometric function, the winding of d around the merinthon is unaffected by a continuous 
gauge transformation. The location where a = 7 r /2 or a = 0 may, however, be shifted by a small 
amount. The merinthons , of course, are not d-function peaks, but extended over a small area. 
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4-5. Wilson Loop Dependence of the Merinthon Solution. 

Each merinthon is constructed from a particular configuration of the 9 field, and each 9 field is 
given by the eigenvectors of a particular set of nested Wilson Loops. This means that in addition 
to the gauge field, the particular topological objects we find will depend on the choice of which 
set of Wilson Loops is being studied. So if we choose one particular set of Wilson Loops, we will 
uncover one distribution of merinthons, while if we choose a different set of Wilson Loops, we may 
well uncover a different distribution of the objects. An obvious question is whether this eliminates 
the usefulness of the merinthons as an explanation for confinement. After all, an instanton, Dirac 
magnetic monopole, or vortex are dependent only on the gauge field, and so can be measured 
without making any arbitrary choices concerning the 9 field. The global topological charge of the 
field is an observable which is independent of any choices that we make. If this can be reduced to 
a set of local topological structures, then the natural thing is to say that these structures should 
also be independent of any arbitrary decisions we should make when parametrising the gauge field 
(such as our choice of Wilson Loop). 

Obviously this issue would not be relevant if the distribution of merinthons was (to a good 
approximation) independent of our choice of Wilson Loop. However, our numerical tests indicate 
that this is not the case: we see no significant correlation between the topological field strength 
when we use one set of Wilson Loops to construct 9 compared to when we use a different set of 
Wilson Loops. 

The Wilson loop represents the propagation of a static quark/anti-quark pair through time. 
The choice of Wilson Loop is equivalent to choosing where the quark/anti-quark pair under study 
is located. So the question posed in the first paragraph of this subsection is asking whether the 
topological objects which explain the confinement of quarks (if quark confinement is explained 
by topological objects) should manifest themselves independently of the presence of those quarks. 
There is no obvious reason why that should be the case. Discussion of confinement makes no sense 
unless there are quarks present to confine; the study of confinement in general can be reduced to 
the study of why each particular pair of quarks is confined. As the explanation of the confinement 
of a particular pair of quarks, where the position of the quarks determines the Wilson Loop 
which then fixes the 9 we require for the decomposition and uniquely determines the topological 
structures contained within 0, the topological objects under study here work as well as those which 
are present in the gauge field independently of any choice of Wilson Loop. There is no obvious 
reason why the key quantity that determines the strength of the confining potential (which is an 
observable independent of which quarks are being studied), the density (after averaging over the 
gauge fields) of merinthon - anti-merinthon pairs, should be dependent on which Wilson Loop is 
chosen. Thus we feel that the arbitrariness of the merinthon fields does not weaken this model’s 
use as an explanation of quark confinement. 

5. String breaking 

So far, we have not discussed the effects of the discontinuity in 9 when the Wilson Loop has 
degenerate eigenvalues. Although the main focus of our work is elsewhere, some discussion of this 
topic is required for completeness; however we will do no more than present a simple result leaving 
the questions arising for the physical consequences of this result unanswered. We will find that 
if the Wilson Loop has degenerate eigenvalues, then the confining potential will weaken at larger 
distances. This is similar to the effect known in dynamical simulations (including quark loops) 
known as string breaking. It is usually expressed in terms of the shielding of the potential by 
the quark loops. Once the quarks become separated far enough, then the binding energy becomes 
sufficiently large that it is advantageous to pull a quark anti-quark pair out of the vacuum, and have 
two shorter ‘strings’ plus an extra pair of quarks rather than one really long ‘string’. Obviously, 
our numerical results are in a pure gauge theory, and we do not expect to observe string breaking. 
If our model is correct, then there must be something suppressing the occurrence of degenerate 
eigenvalues of the Wilson Loop in pure gauge theory, with this suppression weakened in the 
dynamical theory. However, this statement remains just a hypothesis, and a dedicated study, well 
beyond the scope of this work, is required to either confirm or deny it. 
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We expect that the same model of confinement will be valid in both quenched QCD (Pure 
Yang-Mills theory) and full QCD. This means that if confinement is explained by some partic¬ 
ular mechanism just involving topological objects in the gauge field, then there should be some 
explanation in full QCD in terms of those objects or the theory surrounding them why the string 
is broken at large distances. But given that the explanation itself does not depend on the pres¬ 
ence of dynamical quarks - only the weighting given to various gauge field configurations - then 
there should be some allowance for string breaking within the model, and this will apply for both 
quenched and full QCD. String breaking would then be activated by certain gauge fields which 
are suppressed in quenched QCD but present in full QCD. That our model has an allowance for 
string breaking even though we only consider the pure Yang Mills theory is thus not surprising. 

We wish to evaluate the change of the individual gauge links of the Wilson line as we gradually 
increase the spatial extent of the Wilson Loop. Again, we consider a T x i? rectangular Wilson 
Loop, of length L = 2{R -\- T) and T > R. We parametrise the position along the Loop using 
0 < s < L and the spatial extent of the loop with r = R/2. We again neglect the effects of the 
corners of the Loop. 

The building block of the Wilson Loop is 



(90) 


writing 0 = (neglecting the U(l) phase), we can use equation (IA.30I) to write 





n>0 


(91) 


where we use the notation 


[X"Y] = [A, [A,...[A,Y]...]]. 


(92) 


n terms 


a and (j) were chosen so that is Hermitian, traceless and diagonal. 


We now consider the evolution of u with r, using the machinery developed in section [3] and 
equation dSil). We write 



(93) 


with 



(94) 



Using equation (IA.30|) . we can again write this in the form 



(96) 


A2 =5rdrA^ -f ^ ^6r [(-m(/))”a^A] 


(97) 


n>0 
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We can evaluate the last term in equation (IMl) using equation (IA.33I) . which gives 

+ E = cos(2^)(Ai - irtrrAi) + ^[r, Ai] + irtr(rAi). (98) 

n>0 

Given that Ai oc (i.e. one of the diagonal Gell-Mann operators; in SU(2)), we can see that 
tiTAi = 0 and tr(A-l [F, Ai]) = 0. This means that from these two terms only the ^tr A^ cos(2/3)Ai 
part contributes to In most situations, /3 is infinitesimal, and this just leads to = 

-f ... where the ... comes from the remaining terms in equation (1961) and is of 0{6r). 
is Ug plus a small amount, and therefore its expected absolute value increases (as we have seen in 
the previous sections) with r once it is clear of the Wilson Loop eigenvalue degeneracy at r = 0. 
However, where there is a degeneracy in the eigenvalues of the Wilson Loop, /3 will be large, and 
at some point close to the degeneracy cos 2/3 will cross zero. While the first terms in equation (1961) 
will no longer all be infinitesimal, will lose its direct dependence on Ug’L Effectively, the 

value of the Wilson Loop is reset and loses all knowledge of what occurred at i? < 2r. 

Let us suppose that for a particular configuration, these eigenvalue degeneracies in the Wilson 
Loop occur at distances R = 0, Rq, Ri,.... In this case, once we are clear of i? = 0 we may expect 
u (on average) to rise linearly with R until we reach Rq, where it will be reset to some value u. 
Subsequently, it will rise linearly again with R — Rq until R = Ri and it is reset again before rising 
with R — Ri. The location of i?o, Ri etc. will vary from configuration to configuration, which 
means that when we average over configurations, we can consider some mean separation between 
those Wilson Loops with degenerate eigenvalues {Rq). The pattern we may expect is, apart from 
some non-linear behaviour at small R as the loop escapes the degeneracy at i? = 0 the potential 
will increase linearly while R <C {Rq). For R ^ {Rq) we may expect the potential to be flat, 
since on each configuration it will be somewhere between {u) and {u) + p(i?o), and on the average 
over configurations this will become a constant as long as the distribution oi Rq,Ri, ... is not so 
sharply peaked that they occur in the same place on each configuration. There will obviously be 
some intermediate region between these two regimes where we transit from the linear to constant 
behaviour. 

This, of course, precisely as we expect in dynamical QCD: at large distances, the quark loops 
screen the static quark potential leading to a breaking of the string. This has been observed 
numerous times in lattice QCD, for example in [^ . although in our lattice study which excludes 
quark loops we do not see it. 


6. Numerical results 


We generated 16^32 and 20^40 quenched (i.e. pure Yang Mills) lattice QCD (SU(3)) configu- 

3l|, [s^ using a Hybrid 


rations with a Tadpole Improved Luscher-Weisz gauge action [iiliiS 

Monte Carlo routine [s^ (see table |3]) . The lattice spacing was measured using the string tension 
P ^ (420MeV)^. We fixed to the Landau gauge and applied ten steps of improved stout smear¬ 
ing [3J, [33 with parameters p — 0.015 and e = 0. 6 and U were extracted from the gauge field 
numerically by solving equations dzD, (® and (ED, as outlined in [Appendix D| We constructed 
9 from planar Wilson Loops in the xt plane. F^t = is therefore the component of the field 


strength tensor responsible for confinement. All our error analysis was performed using the boot¬ 
strap method (except for the parameters obtained from fits which estimated the errors from the 

distribution). 

The main difficulty we have faced while investigating this proposed mechanism of confinement is 
the non-gauge invariance of 9. We need to work in a gauge where the gauge field is continuous, while 
on the lattice there is no such gauge. However, to completely establish the proposed mechanism of 
confinement, we would need to calculate the variables a and c used to parametrise 9. Extracting 
the components a and c from 9 is straightforward, and we presented some results for the winding 
of c around the peaks in [1^, with some encouraging results. They are, however, gauge dependent, 
and there is therefore no unambiguous definition of these quantities. While the winding number 
is protected from continuous gauge transformations in the continuum, on the lattice it is difficult 
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Lattice size 

L (fm) 

/3 

a (fm) 

# 

16®32 

2.30 

8.0 

0.144(1) 

91 

16®32 

1.84 

8.3 

0.115(1) 

91 

16®32 

1.58 

8.52 

0.099(1) 

82 

20®40 

2.30 

8.3 

0.115(1) 

20 


Table 3: Parameters for our simulations: the lattice size, measured in units of the lattice spacing, the spatial 
extent of the lattice, L, in physical units, the inverse gauge coupling the lattice spacing a, and the number of 
configurations in each ensemble 7 ^. 


even to ensure that you are in the right gauge. Therefore for this paper, we have concentrated on 
gauge invariant quantities such as the static potential and the restricted field strength. 

Our aim is to show that the string tension is reproduced by that of the restricted field, and the 
restricted field dominated by the contribution from 9, or the merinthon term in the field strength 
(constructed from the potential tiXW^dfiO). We should show that the restricted field strength Fxt 
is dominated by points (which would appear on the lattice as peaks a single lattice spacing across), 
while the other components of F contain structures extended in the directions specified in table 

HI 

Measuring the string tension associated with the restricted potential is straight-forward, since 
it is gauge invariant, so it can be extracted from the Wilson Loop using standard methods. Equally, 
the restricted field strength can be measured using the plaquette definition 

_ elu^{x,t)Ut{x + a,t)U^^{x,t + a)U}{x,t)ex, (99) 


and as this is gauge invariant it can be measured without any ambiguities or difficulties. Extracting 
Fj{^, from however, can only be done modulo 2tt for F^ and 47r/-\/3 for F®. For smooth 

gauge fields, where F is small (as is the case in most lattice applications), this is not a concern. 
We, however, as interested in those cases where F is not small, and we shall find that in some 
applications jF-^ j may be larger than tt. For example, if F® ^ 2tt it is impossible to distinguish 
this large value from a fluctuation around zero. This is not of significance for the physics, since 
these shifts in F will not affect the value of the Wilson Loop, but it may be of importance in some 
of the following as we attempt to identify peaks in the field strength. 

However, constructing the merinthon contribution to the restricted field strength is more chal¬ 
lenging because there is no clean gauge-invariant way to observe this on the lattice. In the con¬ 
tinuum, the observable is hmQ_>o {91.6x+ap. — 1) lattice spacing 

a. Since this is not gauge invariant, a gauge transformation which would be discontinuous in the 
continuum would lead to additional discontinuities appearing in the observable, or the removal 
of the discontinuities we want to observe. Fixing the gauge would not help us, because we could 
well be fixing to a discontinuous gauge. We solved this problem by introducing a gauge invariant 
quantity which measures the merinthon contribution to the potential. 

We introduce a field U, which is the gauge field U subjected to a large number of coarse stout 
smearing sweeps. A small amount of smearing is useful to remove discontinuities on the order of 
the lattice spacing; a large amount is usually seen as dangerous as it removes the physical features 
of the gauge field; however we can use this to our advantage. Each Stout smearing sweep replaces 
Ux.fj, —>■ = e'^^^Ux^fj. where Q is a Hermitian operator constructed from closed loops of gauge 

links starting and finishing at x containing the gauge link Ux. n, (t he original smearing algorithm 341 


used plaquettes; we also used 2x1 rectangles following [S^)- After a gauge transformation 


Ux 


^xUx,^,JXx+p, the smeared link transforms in the same way, t/j, ^ > A^ 


TT' A 


The 


second effect of stout smearing is to smooth the gauge field, and reduce the value of the plaquette to 
zero and thus the potential to a constant, and in practice the potential becomes zero. Eventually, 
after repeating the smearing a very large number of times, the over-smeared link U becomes, in 
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effect, a gauge transformation of something close to the identity operator!^ Thus 9l.Ux,^0x+ii is 
both gauge invariant and in a gauge where U is continuous would measure the observable we 

need. We therefore use this observable to measure the merinthon portion of the field strength. In 
effect, we are replicating the calculation for the restricted field, but with the gauge field replaced 
by something close to zero, just leaving the merinthon contribution to the restricted field. To 
extract the Abelian component of OyJx^^Ox+jj, we perform another Abelian decomposition and 

extract the restricted field Ux,ti which satisfies \Ux,fnX^\ = 0 and tr(A'^(A — Al)) = 0 with 

X = Ol.IJx,ij,dx+fi{Ux,iJ.)~^■ We can then compare the field strength from this restricted field, which 
represents the merinthon portion of the restricted field strength, with that of the restricted field 
U. Our expectation is that the observables calculated from the merinthon held and the restricted 
held should be similar: the string tensions should be in agreement, and the held strength should 
contain a similar pattern of peaks. We denote as the held strength tensor constructed using 

(j. 

6.1. String tension - Fixed Wilson Loop 

Firstly, in hgure[2]and table IH we extract the string tension, p, for the original gauge held U, 
the restricted held U and the merinthon component of the restricted held M. We have calculated 
this in two ways: hrstly by calculating the expectation value of the R x T Wilson Loop for one 
of the helds, and htting its logarithm hrstly to VT + a + b/T and then V to pR + c + d/R for 
additional ht coefficients a, b, c and d; and secondly by htting the whole data set according to 
pRT + aR + bT+c+d/R + e/T+ f/(RT) + gR/T+hT/R. The results from these two approaches 
were in good agreement. We tabulate the results from the combined RT ht, while plot the results 
of the single hts. The errors on the expectation value were calculated using the bootstrap method; 
we calculated the errors on the htted parameters by measuring the surfaces of constant In 
principle we should recalculate a new 9 held for each Wilson Loop; however the numerical cost 
to do so was prohibitive for an initial study, and we hrstly show results obtained by using the 
same 9 held calculated from one set of Wilson Loops in the xt plane for our entire simulation [l^ . 
We discuss this simplihcation a little more in the next subsection. We average over all Wilson 
Loops for the restricted held, and not just those with the correct 9 held. Were we to recalculate 
a diherent 9 for each Wilson Loop, the string tension extracted from the U and U helds would 
be identical. Because we have not done this, we expect the two quantities to differ by a small 
amount. 

We calculate U after 100, 300, 500, 600, 800 and 1000 sweeps of stout smearing with parameters 
e = 0, p = 0.015 following the notation of [35l| . and extract the merinthon portion of string tension 
by constructing Wilson Loops from the Abelian decomposition of 9^U9. We also show the string 
tension for {7, and can conhrm that it is much smaller than that of the original gauge held and 
decreases as we increase the level of smearing. The observable we are using to measure the 
merinthon portion of the string tension is unaffected by the smearing, suggesting that we have 
indeed measured the contribution from 9dfj,9^ rather than any remnant of the gauge held remaining 
after the smearing. We see a good agreement between the merinthon and restricted held string 
tensions, within 10% on the smaller lattices (although the difference is larger on the larger lattice 
where we have less statistics), suggesting that it is indeed structures within the 9 held that lead 
to conhnement. There is no obvious change in this pattern as we change the lattice spacing. 

6.2. String Tension - Varied Wilson Loop. 

We generated the results presented in the previous section as a quick and dirty initial study, 
largely to demonstrate to ourselves and the wider community that this method was worth pursuing, 
and we hrst published these results in [1^ (these are the 10 smearing sweeps results of table [5]). 


®It has been shown that a large number of infinitesimal smearing sweeps - not quite what we are doing here, 
but similar - is equivalent to a flow in the configuration space towards an ensemble generated with infinite gauge 
coupling [s^l . 
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0 = 8.52 



Figure 2: The string tension extrapolated to infinite time for the original gauge field f/, the restricted gauge field 
U^ the over-smeared field U and the merinthon field M for the (3 = 8.52 ensemble 


0 

8.0 

8.3 

8.52 

8.3L 

u 

0.094(2) 

0.0590(8) 

0.0442(6) 

0.057(2) 

u 

0.116(4) 

0.095(2) 

0.077(1) 

0.099(1) 

Uwo 

0.0828(4) 

0.0545(3) 

0.0433(7) 

0.0594(8) 

ATioo 

0.129(4) 

0.090(3) 

0.075(3) 

0 .100(1) 

U300 

0.0460(3) 

0.0301(3) 

0.0239(4) 

0.0316(2) 

M300 

0.122(4) 

0.086(3) 

0.072(3) 

0 .102(1) 

U500 

0.0316(2) 

0.0216(2) 

0.0174(3) 

0.0218(1) 

M500 

0.124(5) 

0.088(3) 

0.072(3) 

0.104(1) 

Uqoo 

0.0273(2) 

0.0185(1) 

0.0149(2) 

0.0179(1) 

Mqoo 

0.103(10) 

0.087(5) 

0.076(3) 

0.099(1) 

Usoo 

0.0213(1) 

0.0147(1) 

0.0124(2) 

0.0144(1) 

Msoo 

0.104(8) 

0.087(5) 

0.076(3) 

0.099(1) 

Uiooo 

0.0174(1) 

0 .0122(1) 

0.0106(2) 

0 .0121(1) 

Miqoo 

0.105(9) 

0.087(5) 

0.077(3) 

0.098(2) 


Table 4: The string tension extrapolated to infinite time across all our ensembles. 8.3L refers to the 20^ x 40 
ensemble. 
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Since the initial drafts of this paper, we have commenced a more accurate study removing the 
simplification of the previous section, this time recomputing the Abelian decomposition for each 
Wilson Loop under consideration. This obviously maintains the identity between the confining 
potential of the full and restricted gauge fields. Our results so far are only tentative, and far 
from being finalised. We present early results here mainly to update our initial presentation of 
this study in [s^- In those proceedings, we found that although the restricted and full QCD 
potentials were in perfect agreement, the topological contribution to the string tension seemed to 
be relatively small. 

If this result is genuine, there are two possible explanations: firstly, that the Maxwell term in 
the restricted potential, contrary to expectations, also has a large contribution to confinement. 
The second possibility is that our means of extracting the topological part of the string tension 
is flawed, and what we were measuring was not the topological part of the string tension but 
something else. Since our initial publication, however, we have expanded our study, and have 
noted another issue. We apply a small number of smearing sweeps to the gauge field initially, 
before applying the Abelian decomposition, to remove various unphysical fluctuations of the order 
of the inverse lattice spacing. This is necessary to allow us to obtain a clear result. How many 
smearing sweeps should be used is difficult to assess, and for a full presentation of the results 
we should measure the string tension at different levels of smearing to gain some idea of the 
systematic error caused by this arbitrary choice. As we did so, we observed that the measured 
topological contribution to the smearing depended strongly on the number of initial smearing 
sweeps. Furthermore (unlike in the results of the previous section), we also found that at a 
larger number of initial smearing sweeps, there was a strong lattice spacing dependence on the 
ratio between the restricted and topological potentials: as the lattice spacing grew smaller, the 
restricted potential became increasingly dominated by the topological part. This means that our 
full analysis is considerably more difficult to do accurately than we had first envisaged, and, given 
our limited computing resources, will take longer to complete. 

Some preliminary results (which are not finalised, and could change when we make our final 
publication) are shown table [3 As we vary the number of smearing sweeps, the string tension of 
the original gauge field first increases as the unphysical dislocations are removed, and then slowly 
decreases as the smearing gradually distorts the gauge field. This obviously begs the question of 
which value we should use, and our intention has been to use the number of smearing sweeps where 
the string tension is at its maximum (at the present time, we have kept the smearing parameters 
fixed at e = 0, p = 0.015 for this study). So far, we have only conducted this analysis for the 
16^ X 32 configurations. Table [5] lists the /3 value of the configuration, the number of smearing 
sweeps where the maximum value of the string tension was reached (out of those which we have 
sampled so far), and the measured string tension for the QCD (and restricted) gauge fields, and 
the proportion of that string tension which we have measured as coming from the topological part 
of the restricted potential (after 2500 smearing sweeps). 

It can be seen that firstly the maximum string tensions for the /38.3 and /38.52 configurations 
are higher than the results given above (they were calculated following a different amount of 
initial smearing). Secondly, the proportion of the string tension attributable to the topological 
part increases as we decrease the lattice spacing. This pattern is not merely a matter of the 
results being taken after a different number of smearing steps, but is apparent at each value of the 
initial smearing, although it is more pronounced as we increase the smearing. Our analysis is less 
developed on the larger lattice, but the results we have suggest that there is no significant volume 
dependence. 

This strong dependence on the lattice spacing is awkward because it necessitates an extrapo¬ 
lation to the continuum limit. This is challenging with the data we have gathered so far, firstly 
because we do not know of a good formula to guide the extrapolation, and secondly because we 
are still on a fairly coarse lattice. We need to finalise our results on these small lattices, and then 
expand our data set so that we have access to hner lattice spacings before it is useful to attempt 
to perform such an extrapolation; we will have to leave the completion of this study to a future 
publication. 

However, it is not unreasonable to suggest on the basis of this data that the string tension will 
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p 

8.0 

8.3 

8.52 

Smearing Sweeps 

P 

Proportion 

10 

0.0976(36) 

0.254(14) 

10 

0.0603(19) 

0.277(14) 

10 

0.0416(17) 

0.334(19) 

Smearing Sweeps 

P 

Proportion 

16 

0.0954(40) 

0.247(17) 

16 

0.0740(17) 

0.332(15) 

16 

0.0489(14) 

0.462(27) 

Smearing Sweeps 

P 

Proportion 

22 

0.0944(26) 

0.334(15) 

22 

0.0.0731(14) 

0.400(14) 

22 

0.0501(13) 

0.526(27) 


Table 5: The number of smearing sweeps where the string tension of the full QCD field was at its maximum, the 
string tension p (calculated from the U field) and the ratio of the topological part of the string tension to the full 
QCD string tension. Only the statistical errors are recorded. We list results at those numbers of smearing sweeps 
where, out of the measurements we have taken so far, the string tension for the /98.0, /38.3 and /38.52 ensembles 
respectively are maximised. 


be dominated by the topological part of the restricted potential (as we have measured it) in the 
continuum limit. This still needs to be confirmed by a complete and more careful analysis, which 
we are currently undertaking. 

6.3. Presence of peaks in Fxt 

The next question is whether the component of the gauge invariant restricted field strength 
is dominated by point-like objects as expected in the theoretical model, and we plot the distribution 
oiF^, in figure [3] on two neighbouring slices of the lattice, using a contour plot to display the data. 
The corresponding plot for has a similar set of structures. It can be seen that is indeed 
dominated by these objects a single lattice spacing across (or, on occasion, two lattice spacings - 
see also figure (HI). There is no relationship between the peaks on neighbouring slices of the lattice, 
suggesting that these are indeed points rather than strings or surfaces. Are these peaks from the 
merinthon part of the restricted potential? The background shading of figure [3] gives a similar plot 
for the field strength extracted from 17, and the is a strong correspondence between the location 
of the peaks in the merinthon field and the restricted field, although on occasion one of the peaks 
may be shifted by a lattice spacing (which is the limit of the resolution of the lattice, so we can 
expect errors of up to a half lattice spacing on each of the data sets), and a few objects in the 
merinthon field are not visible in the restricted field. Nonetheless, the similarity between the two 
sets of data indicates that these peaks are caused by the magnetic portion of the ensemble as 
predicted by the theory. Although figure [3] shows data from just one slice of one configuration for 
one ensemble, the pattern of peaks is similar across all our ensembles. 

In figurelHand tablelH we investigate whether these peaks are responsible for the string tension, 
i.e. if excluding the peaks would reduce or eliminate the confining potential. Usually, we measure 
the expectation value of the Wilson Loop by averaging over every planar loop in the configuration. 
However, it is a straight-forward exercise to restrict the average to loops in the xt plane and either 
only include or exclude those Wilson Loops which contain one of the peaks from the average. We 
show in figure |4] the result of excluding the Wilson Loops which contain peaks. We introduce a 
cut-off C, and do not include any Wilson Loops where |F| > C for one of the lattice sites contained 
within the loop. If the peaks are responsible for confinement, we would expect to see the string 
tension disappear as we exclude more of the peaks. If they are merely incidental, and confinement 
is dominated by some other mechanism than that proposed here, we might expect to see little 
difference in the string tension as we change the cut-off. The picture is, however, consistent with 
the idea that these peaks are responsible for confinement. 
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Figure 3: A comparison between the peaks in (contours) against the merinthon field strength extracted U 
(shaded background). In this presentation the negative and positive peaks cannot be distinguished. The top plot 
is taken for slice of the lattice Y = 0, Z = 10 while the bottom plot uses the slice at Y = 0, Z = 11 on the same 
(3 = 8.52 configuration. Due to the limited resolution of the lattice, the extrapolated contour lines and the shading 
have an error of up to one lattice spacing. 


P = 8.52 Peaks in Fxt excluded, T = oo 



Figure 4: The U string tension, p, excluding Wilson loops containing peaks of height \Fxt\ > C from the average, 
on the ^8.52 lattice. 
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c 

2.55 

2.30 

2.05 

1.80 

1.55 

p = 8.0 

0.1178(9) 

0.1201(9) 

0.1213(9) 

0.1225(8) 

0.1199(8) 

^ = 8.3 

0.0943(8) 

0.0955(7) 

0.0961(8) 

0.0971(8) 

0.0957(8) 

^ = 8.52 

0.0796(7) 

0.0809(7) 

0.0820(6) 

0.0821(6) 

0.0820(6) 

/3 = 8.3F 

0.0975(8) 

0.0965(8) 

0.0954(8) 

0.0964(8) 

0.0963(8) 

C 

1.30 

1.05 

0.80 

0.55 

0.30 

P = 8.0 

0.1144(8) 

0.1049(7) 

0.1002(7) 

0.0925(7) 

0.04998(7) 

P = 8.3 

0.0899(8) 

0.0866(7) 

0.0809(6) 

0.0741(6) 

0.0382(7) 

P = 8.52 

0.0791(6) 

0.0747(6) 

0.0697(5) 

0.0665(5) 

0.0472(6) 

/3 = 8.3F 

0.0902(7) 

0.0843(7) 

0.0770(7) 

0.0706(5) 

- 


Table 6: The U string tension, p , excluding Wilson loops containing peaks of height \Fxt\ > C from the average. 
The /38.3L data refers to the largest lattice; there was insufficient data to get a reliable estimate of the string tension 
at the lowest cut-off on this lattice. 


6.J^. Electromagnetic field tensor 

We now look at the full strength tensor to examine how closely it resembles our expectations 
from section 14^31 although for conciseness we concentrate only on (the results for and 
are indistinguishable). First of all, we plot histograms of each component of the field strength 
in figure [5l showing what proportion of lattice sites had a particular value of the electric or 
magnetic field. It can be seen that the distribution for Ey, By, E^ and are all similar, while the 
distribution for Ex and Bx differ. Most of the plots show similar features: a large peak around 
H^i, = 0, meaning that most of the points have no field strength beyond a small fluctuation around 
zero, and a plateau for larger field strength up to the maximum value of meaning that on 
a minority of lattice sites there is a large contribution to the field strength, where the number of 
sites with that particular value of the field strength seems to be independent of the field strength. 
This is as we expected. We also see that the proportion of the points with small values of the 
field strength increases as we decrease the lattice spacing, which suggests that these objects have 
a fixed density per physical volume but the size of the objects decreases as the lattice spacing 
tends to zero, which is consistent with the idea that these are thin points or string-like objects. 
Bx and Ex are different; there are twice as many lattice sites with a large value of Bx compared 
to the fields in the y and z direction, and very few lattice sites with a large value of Ex- There 
also seems to be little volume dependence. 

The picture for all the fields is consistent with having some (perhaps Gaussian distributed) 
fluctuations around zero, and then some additional objects at large field strength. For Ex, the 
number of the objects diminishes as the field strength increases, while for the remaining fields the 
proportion of lattice sites with a particular field strength seems to be roughly independent of the 
field strength, at least for those field strengths we can measure. This means that there may be 
peaks in the field with a field strength ^ 27r for F® or ^ 47r-\/3 for F® which are impossible for 
us to distinguish from zero. It should therefore be borne in mind that not every object in these 
fields will be visible to us. 

This picture can be studied in more detail by plotting slices of the electromagnetic fields: figure 
[6] shows this in the xt plane; figure [3 in the yt plane, [8] in the zt plane. 

These objects are not just isolated magnetic monopoles: there are peaks in the electric field, 
and the magnetic fields are not spherically symmetric. They may, of course, be some magnetic 
monopoles and some other structures, or possibly a signature of a monopole/anti-monopole con¬ 
densate. The results are consistent with our expectations given in table [U and the patterns seen 
on this configuration and slice of the lattice are duplicated across all our ensembles. The electric 
field is clearest (although not completely unambiguous): Ex is a point, while Ey and E^ form 
strings parallel to the a;-axis. The magnetic field is a little less clear, and expected string like 
behaviour of By and Bz is less obvious, although we have seen clear examples of t-strings. Bx is, 
as expected, extended along both the x and t axes, though not the z or y axis. In particular, other 


33 








Distribution Distribution 



Field strength 



Field strength 


o 


^3 


O 


X) 


P 


0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 
0 

-2-101234 
Field strength 

0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 
0 

-2-101234 
Field strength 




Figure 5: Histogram of the various components of the magnetic field strength for /9 = 8.0 (top left), jS = 8.3 (top 
right) and j3 = 8.52 (bottom left) on the 16^32 lattices, and the /3 = 8.3 20^40 lattice (bottom right). 
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Figure 6: The electromagnetic fields on the slice with ^ = 5, 2 = 5 on one (3 = 8.52 configuration. The dotted 
contours indicate a positive field strength, the solid contours indicate a negative field strength. The top plots show 
the fields in the x direction, the middle plots the y direction, and the bottom plots the 2 : direction, while the left 
column shows the electric field and the right column the magnetic field. 
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Figure 7: The electromagnetic fields on the slice with a; = 5, 2 : = 5 on one (3 = 8.52 configuration. The dotted 
contours indicate a positive field strength, the solid contours indicate a negative field strength. The top plots show 
the fields in the x direction, the middle plots the y direction, and the bottom plots the 2 direction, while the left 
column shows the electric field and the right column the magnetic field. 
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Figure 8: The electromagnetic fields on the slice with a: = 5, y = 5 on one /3 = 8.52 configuration. The dotted 
contours indicate a positive field strength, the solid contours indicate a negative field strength. The top plots show 
the fields in the x direction, the middle plots the y direction, and the bottom plots the 2 : direction, while the left 
column shows the electric field and the right column the magnetic field. 
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orientations of the fields with strings in directions perpendicular to the xt plane are inconsistent 
with this data. 


6.5. Cluster analysis 

To investigate the restricted electromagnetic fields more thoroughly, we employ a simple strat¬ 
egy to determine the size and shape of these structures. We first of all find clusters, connected 
regions of sign coherent high electric or magnetic field strength. We perform the analysis for each 
component of the field strength tensor separately. While we have performed this analysis for 


^ /iI/5 


and individually, there is no significant difference in the results, so here we just 


show the results for 

We dehne a cluster as all connected sites with a value of \F^v\ > 1 and identical sign of F^^ (i.e. 
except for clusters containing just one lattice site, each site in the cluster is the nearest neighbour 
of at least one more site in the cluster and all its nearest neighbours which satisfy the bound 
\F^y \ > 1 are within the same cluster and each cluster is internally connected so it is possible to 
reach any site within the cluster from any other by crossing to a neighbouring site while remaining 
within the same cluster, while different clusters are not internally connected with each other). We 
have also studied (but do not show here) results with numerous values of this cut-off from 0.8 to 
2 .0, but we do not see any important difference as a result of varying this cut-off. 

For all the fields except Ex, the vast majority of the lattice sites where the restricted field 
satished our condition were contained within a large cluster which spanned the entire lattice. This 
is in certain respects similar to the structures found in the full field strength tensor 38, 3^. The 
theoretical picture we were testing described a number of strings in either the x ov t direction. 
This large structure may be formed if the density of these strings is sufficiently large that every 
string has at least one neighbour on an adjacent lattice site in the y oi z directions. It may, 
however, relate to a different picture of the vacuum. There were a few clusters not attached to 
this large structure; and most of our analysis is dominated by results from these smaller structures, 
as our averaging with respect to the number of structures weights against the smaller number of 
structures with a large number of lattice sites. We also present a few results which specifically 
target these larger structures. Here we see the expected string behaviour, with a large number of 
objects at just a single point. We emphasise that Ex (or Fxt) differed, in that the field strength 
did not contain this large structure, but only single lattice site structures. To attempt to analyse 
the shape of this large cluster, we have also performed a similar analysis, but where the clusters 
are constructed from connected sites with a two dimensional slice of the lattice. 

When we show our results, we will label those plots where we have included all clusters as 
unfiltered, those which only average over clusters with at least four lattice sites as filtered, and 
those which average over clusers with at least two hundred lattice sites as Large Clusters. 

The theoretical expectation was that if Ex is dominated only by points of large field, then the 
remaining fields should be dominated by objects which are either points or strings in the x or 
t direction. The only expected difference in the results for the Ey, Ez, Bz and By fields is the 
orientation of the string. A point will show itself as a cluster containing just a few lattice sites 
(ideally only one lattice site, but due to lattice artefacts it may be smeared across two or three 
sites). A string would display itself as a line of flux, i.e. each site within a cluster would have only 
two nearest neighbours within a cluster if it is in the centre of the string and one if it is at the 
end of the string, and it will only be extended in one direction. In practice, the situation will be 
messier than this. Neighbouring strings may be counted as the same cluster, and there may be 
ambiguities caused by the poor resolution of the lattice leading to a thickness of the string being 
more than one lattice spacing in places, two different strings may overlap (and the expectation 
is only accurate if the functions G and R' have no strong angular dependence). Nonetheless, we 
would hope to see some signal of these strings when considering an average size and shape of the 
cluster. 

We first plot the distribution of the number of lattice sites in each cluster in figure [9] The 
picture is similar across all our ensembles, so we have only plotted the data for the /3 = 8.52 
ensemble. It can be seen that all the components of the field strength contain a large number of 
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Figure 9: The distribution of the number of the lattice sites in each cluster for the /3 = 8.52 ensemble. 

points with only one or two sites contributing to the cluster. However, the field only contains 
these objects, while the remaining fields also contain a substantial number of structures extended 
across several lattice sites. This means that if the picture outlined in the earlier sections is correct, 
then the Ey, E^, By and B^ fields can only contain points and the strings described in table 
HJ We also note that the data for the Ey, E^, By and Bz fields are indistinguishable, in line with 
our expectations. 

We next study the number of nearest neighbours within a cluster in figures [TUI ttH and dH 
The first of these plots shows the data across all the clusters, while the second filters only those 
clusters containing four or more lattice sites, so eliminating the points. We have studied the effects 
of changing this filter from all clusters with at least three sites to all clusters with at least six 
lattice sites, and there are no substantial changes in the results. The third plot focusses on the 
very large clusters. A point should have no nearest neighbours (or possibly on occasion one), a 
string one or two neighbours, a two dimensional surface one, two three or four neighbours (with 
the proportion depending on the size and shape of the surface), and higher dimensional objects 
would have a larger number of nearest neighbours within the cluster. The data for E^ is again 
unambiguous: these are points. For the remaining fields, it is less clear. Focussing on figure [TTl 
we see that firstly there is no substantial difference between the results for the different fields: the 
same structures appear in the fields with the same frequency. The majority of the lattice sites 
within the larger clusters had two or three neighbours, with a small number having four. We can 
rule out that these are extended in more than two dimensions. This may indicate strings where 
two strings are on neighbouring lattice sites and thus counted in the same cluster (and giving 
objects with three or more nearest neighbours), or it may indicate that there are also some two 
dimensional surfaces with a large electromagnetic field. There is a small but statistically significant 
increase in the proportion of sites with two nearest neighbours as we decrease the lattice spacing, 
and no noticeable difference in the distribution as we increase the volume of the lattice. For the 
very large clusters, which are extended across the whole lattice, we again see that the sites have 
two or three nearest neighbours in the cluster. This is consistent with, but does not demonstrate, 
that these clusters are constructed from densely packed one dimensional objects. 

We predicted that the strings in the B^ field could be in either the x oi t direction, while 
the electric fields Ey and Ez would extend in the x direction only and the magnetic fields By 
and Bz in the time direction only. In figure 1131 we plot the distribution of the spatial extent in 
each direction of each cluster (for example, a cluster that stretched from (t, x, y, z) = (3,4, 5, 7) to 
(6, 9,8,13) would have extents (6 — 3 + 1 = 4,6,4, 7); the plot shows the proportion of clusters with 
a particular extent in each direction). An idealised cc-string (which would not occur in practice) 
would have a large extent in the a:-direction, and extend for only one lattice site in all other 
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Figure 10: The distribution of the number of the number of nearest neighbours within a cluster for each site in 
a cluster for j3 = 8.0 (top left), /9 = 8.3 (top right), /3 = 8.52 (bottom left) 16^32 lattices and the /3 = 8.3 20^40 
lattice (bottom right). 
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Figure 11: The distribution of the number of the number of nearest neighbours within a cluster for each site in 
a cluster for /9 = 8.0 (top left), fi = 8.3 (top right), /3 = 8.52 (bottom left) 16^32 lattices and the /3 = 8.3 20^40 
lattice (bottom right), only including those clusters containing at least four lattice sites. 
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Figure 12: The distribution of the number of the number of nearest neighbours within a cluster for each site in 
a cluster for /3 = 8.0 (top left), j3 = 8.3 (top right), (3 = 8.52 (bottom left) 16^32 lattices and the (3 = 8.3 20^40 
lattice (bottom right), only including those clusters containing at least 200 lattice sites. 
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directions. This is not fully realised in our data: there is a noticeable number of structures which 
are extended for up to four or five lattice sites in the ‘unwanted’ directions. However, for each field, 
the clusters extend considerably more in the directions of the expected strings. We once again 
restrict the results we present to those clusters containing at least four lattice sites, although once 
again we have studied numerous values of this filter without seeing any significant difference in the 
results. The shape of the curve changes as we only include larger clusters (larger extents become 
more likely as we increase the cutoff for the cluster size), although the maximum in the X and 
T directions remains at four lattice spacings. This may be consistent with our expectations, with 
strings overlapping each other, or it may indicate that there are additionally other types of objects 
in the electromagnetic field. 

One question that can be raised is whether that some clusters have a large thickness in the 
unwanted directions and some sites within the clusters challenge whether these large fields are 
arranged in strings one lattice spacing thick. Clearly, in the context of this picture, this can be 
explained by having strings on neighbouring lattice sites, and some degree of the breakdown of 
the naive behaviour should be expected. To better judge whether this is indeed the case, we have 
generated configurations where the field strength is distributed as expected. We have counted 
the number of lattice sites where > 1 and < —1, and generated two different sets of 
configurations of the restricted field strength. In each case, we ensured that the same number of 
lattice sites had the large F^ as on the actual quenched QCD configuration. In the first random 
configuration, these sites were allocated completely random positions around the lattice. In the 
second random configuration, one quarter of the sites were positioned randomly, and the rest were 
arranged into strings in either the x or t direction as appropriate for the field. The principle 
ambiguity is what to use as the distribution for the length of these strings, since the distribution 
given in figure [T^] will contain examples from multiple strings. We examined several Poisson and 
Gaussian distributions, and the data we present here is taken from a Poisson distribution with the 
mean value extracted from the data presented in figure fTSl Only the analogue of the curve in m 
in the direction of the string depends strongly on the distribution of the length of the string. We 
cannot predict the shape of this curve for this simple model. We then simply repeat the cluster 
analysis on these configurations. 

Some sample results are plotted in figures (TT] and [131 

It can be seen that the data for figure [M] is qualitatively similar to the actual QCD data for 
Ex, while there are disagreements when comparing the model and actual data sets for the other 
fields. The second model, with strings, compares well with the numerical data for By, Bz, Ey 
and Ez', in particular the plots for the number of nearest neighbours and the extent of cluster in 
the directions away from the string are close to the QCD data. The extent of the cluster in the 
direction of the string strongly depends on the distribution for the string length, so we do not 
expect a perfect agreement here. In particular, that the ‘thickness’ of the string agrees with the 
data, and seems to be largely model independent, suggests that this cannot be used to rule out the 
predicted string model of flux. However, the broad agreement between our model calculation and 
the data implies that the QCD Ex field is dominated by points and the other fields by a mixture 
of points and strings. 

Finally, we check to see whether the locations of the clusters in the different fields are correlated. 
If our model is correct, then the space-time location of many of the peaks in Ex should be correlated 
with high field strengths in the other fields. To do this, we take each of the peaks in Ex (i.e. each 
point where \Hxt\ > 1) and plot a histogram of the value of the other fields around that site. For 
comparison, we produce a similar plot measuring the field strength around random lattice sites 
rather than the peaks of the Ex held. We would expect the plot for the actual data to have a higher 
concentration of large field strengths than the random data, and a smaller concentration of lower 
field strengths. There is one complication: we measure the field strength at different locations 
in space-time for each field. For example, we measure Fxt{t -\- \,x + \,y, z) = Ex(t,x,y, z) and 
Fyz{t, X + ^,y,z + i) = By{t,x,y, z). This means that there may be a shift of up to one lattice 
spacing in the location of the peak within each of the measured fields. We therefore take the 
maximum value of \By{t, x, y, z)\, \By(t + I, x, y, z)\, \By(t, x,y, z — 1)| and \By{t + l,x,y, z — 1)| 
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Figure 13: The distribution of the spatial extent of the clusters for the /3 = 8.52 ensemble, only including those 
clusters containing at least four lattice sites. 
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Figure 14: The distribution of the spatial extent of the clusters for where the field strength is selected according to 
the distribution given in figure^for the (3 = 8.52 ensemble for the By (top left) and Bx (top right) fields, but with 
the lattice sites randomised; the distribution of the number of nearest neighbours (bottom left), and the filtered 
distribution of the number of nearest neighbours (bottom right) for this random ensemble. 
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Figure 15: The distribution of the spatial extent of the clusters for where the field strength is selected according 
to the distribution given in figure O for the ^ = 8.52 ensemble for the By (top left) and Bx (top right) fields, but 
with the lattice sites randomly distributed with three quarters of the larger field strengths arranged into strings; 
the distribution of the number of nearest neighbours (bottom left), and the filtered distribution of the number of 
nearest neighbours (bottom right) for this random ensemble. 
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(for example). This creates a bias towards larger field strengths for both the actual locations of 
the peaks in and the random data. This is particularly pronounced for where there are 
16 lattice sites which we need to consider, and coupled with the higher density of lattice sites 
with a large value of the Bx field this means that actual data from Bx does not differ greatly 
from the randomised data. We plot the results in figure [16] for the /3 = 8.52 ensemble. There is 
a statistically significant difference between the curves for the Ex peaks and the random data for 
all the fields, with a larger proportion of sites with a large value of the By, Bz Ey and Ez fields 
and a smaller proportion of sites with a small value of these field surrounding the peaks in Ex 
than for the random data. There is, however, still a large number of sites where there does not 
seem to be a great deal of correlation between the peaks in (for example) Ex and Ey] partly this 
can be explained by those sites where in practice Ey ^ Ey/Zir but we measure Ey = 0, but it still 
seems as though not every peak in Ex corresponds to a peak in every one of the other helds. This 
is not inconsistent with equation (l82l) : for example if diG is small then only three of the other 
fields are large when there is a contribution in Ex, and if both diG and d 2 G are small then only 
two of them do. This data is then qualitatively but not conclusively in line with the expectations: 
there is a clear correlation between the fields, however it is not large enough to say that all the 
the topological structures in the different components of are correlated with each other. 

This analysis has concentrated on the smaller, isolated structures. We still have to consider 
the larger structure which extends throughout the lattice, which we have seen from the nearest 
neighbour analysis is either one or two dimensional. This is consistent with our expectations if it is 
constructed from neighbouring strings along the X oi T directions (so we may have an a;-string at 
y = z = t = 0 and extending from a; = 0 to x = 5, and another string at y = 1, z = t = 0 extending 
from X = 5 to X = 10, and so on; if these are closely enough packed then our analysis will place 
them into the same cluster). We investigate this structure by modifying the cluster algorithm. We 
firstly select two directions; and we construct the cluster by only considering nearest neighbours 
in those two directions. This restricts us to isolated objects. By varying the two directions, we 
can see if the strings predominantly extend in particular directions. 

We display the extent of the strings in figures [T71 [T51 and fTOl We have shown the xt, xy and yt 
planes as examples. There was again little difference between our ensembles. Again, we restrict 
ourselves to clusters with at least four lattice sites, so there must be an extent of on average two 
or three lattice sites in at least one direction. We again see that the magnetic fields By and Bz are 
extended far more parallel to the t axis than in the other directions. Ex and Ez along the x-axis, 
while Bx is extended along both the t and x axes. None of the fields show similar features parallel 
to the y and z axis. 

The picture we see is consistently that the restricted Ex field is dominated by points, while 
the Bx field by strings in both directions in the xt plane, and the other electric fields by strings 
in the x direction, while the other magnetic fields by strings in the t direction. These are not 
magnetic monopoles (or not only monopoles): the behaviour of the electric fields and Bx rules 
this out (and note that we presented results for the supposedly ‘magnetic’ trn[dyn,dun\ part of 
the field strength; although the picture for the full field strength is the same). The results for the 
restricted electric and magnetic fields are consistent with the expectation (though, of course, they 
do not prove that the model of confinement we presented in the first part of this paper is correct: 
they may also be consistent with other models of the vacuum). 

7 . Conclusions 

We have investigated whether the confining quark potential in quenched QCD is caused by CDG 
gauge-invariant merinthons. Our main novelty is to construct the CDG colour field = 9X^9^ 
from the eigenvectors of the Wilson loop, which allows us to access the full symmetry group of 
the merinthons and permits a theoretical discussion. From this analysis, we have shown that the 
Wilson Loop for the actual QCD gauge field is identical to the Wilson Loop constructed from the 
restricted Abelian field, and this allows us to express the Wilson Loop in terms of an area integral 
over its surface. We then find that the construction of 9 allows for certain topological defects to 
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Figure 16: The distribution of the value of the field at the location of the peaks in the Ex field and at random 
lattice sites for the /3 = 8.52 ensemble. 
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Figure 17: The distribution of the spatial extent of the clusters in the xt planefor the (3 = 8.0 ensemble, only 
including those clusters containing at least four lattice sites. 
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Figure 18: The distribution of the spatial extent of the clusters in the yt planefor the (3 = 8.0 ensemble, only 
including those clusters containing at least four lattice sites. 
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Figure 19: The distribution of the spatial extent of the clusters in the xy planefor the (3 = 8.0 ensemble, only 
including those clusters containing at least four lattice sites. 
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appear in the restricted field strength. In the component of the restricted field strength in the 
plane of the Wilson Loop, these manifest themselves as point-like objects; in the other restricted 
electric and magnetic fields they manifest themselves as strings parallel to one of the axes of the 
Wilson Loop. The net result is that the restricted field strength is dominated by a large structure 
which is constructed from these one dimensional strings. The full QCD field strength contains the 
restricted field, 


= hAM + FAX] + A]), ( 100 ) 

so it may be expected that these global one dimensional structure may extend to the full QCD 
field strength (albeit with 0(4) symmetry suggesting that F^^,\X] will also contain strings in 
other directions). If these strings are sufficiently dense, this picture is consistent with lattice 
simulations [s^, [s^. In particular, these objects are not magnetic monopoles; at least for this 
particular choice of the colour field. The peaks in the Ex field lead to an area law for the Wilson 
Loop and thus confinement. 

We have also given an argument within the context of this decomposition which might explain 
string breaking at large spatial distances. 

Our numerical results focus on searching for these strings and peaks in the restricted field 
strength, and we have uncovered firm evidence that the electric field in the direction of the Wilson 
Loop, which is responsible for the confining potential, is dominated by objects which are just 
present on a single lattice spacing. The remaining components of the electric and magnetic fields 
contain one dimensional objects which extend in the directions suggested by our theoretical analy¬ 
sis. While this does not prove that our analysis is correct, they are certainly consistent with it; and 
whatever these peaks are, they are responsible for the confining potential. Our numerical analysis 
is consistent with the vacuum related to the restricted field being dominated by a combination 
of magnetic monopoles and the tti topological objects we have labelled merinthons. It is also 
possible that the binding of monopoles and anti-monopoles into pairs (similar to the Cooper pairs 
in superconductivity) might lead to a monopole condensate with a different pattern in the field 
strength than one would expect from a naive picture of isolated monopoles and anti-monopoles, 
and we cannot rule out this possibility. We also show that in an initial test, making an important 
approximation, the string tension can be accounted for by the second, ‘magnetic’, term within the 
restricted field strength. The results of the full calculation (without any approximation) are less 
clear: all we have shown so far is that we will need to make a careful study of the continuum limit. 
The objects responsible for confinement are embedded in the colour field, and only indirectly in 
the gauge field. 

We believe that this evidence is enough to merit further, more detailed, investigation of whether 
confinement can be explained within the context of this particular formulation of the CDG Abelian 
decomposition. 

We have not yet investigated the effects of changing the representation of the group, of incor¬ 
porating quark loops into our lattice study, or studied the de-confinement transition, and hope to 
investigate at least some of these questions in future work. 
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Appendix A. The Abelian Decomposition 

Appendix A.l. Proof of the existence of 9 

In the text, we defined 9 as the field which diagonalises the gauge link 11^ along the Wilson 
Line. We here demonstrate that this field exists, using the following three lemmas (and the well 
known properties of the eigenvectors of normal and unitary matrices, including that a unitary 
matrix can always be diagonalised). 

Lemma Appendix A.l. Suppose that there is a U{Nc) field 9 which diagonalises the gauge 
links Ua which contribute to the closed Wilson Line W[Cs] of length L starting and ending at 
a position s, so 9\U„9„+Sa = for diagonal Gell-Mann matrices Xfi Then 9s contains 

the eigenvectors ofW[Cs\, with eigenvalues , i.e. (given that W is unitary and thus 

normal) W[Ca]9s = 9se'‘^<' Sau^x^ ^ 

Proof W[Cs] is defined as 

W[Cs] = UsUs+Sa ■ ■ ■ Us+L-2SaUs+L-Sa- (A.l) 

and, since 9s = 0s +l as s and s + L refer to the same location in space, 

W[Cs]9s =UsUs+Sa ■ . ■ Us+L-25aUs+L-5a9s 

= UsUsUs+Sa ■ ■ ■ Us+L-2Sa9s-\-L-Sae^^‘+^-^'’^ 

=UsUsUs^&. ■ . 

(A.2) 

where the last line follows by repeatedly evaluating \J9 = and we have noted that the 

diagonal Gell-Mann matrices commute with each other. But, since is diagonal, this is 

just an eigenvalue equation. Hence we have proved the result. □ 

Lemma Appendix A.2. The eigenvalues of the Wilson Line along a closed curve C are inde¬ 
pendent of which location s along the curve is used as the start and end of the line. 

Proof We note that 


W[Cs] =UsUs+5. ■ ■ ■ Us+L-25aUs+L-Sa 
— Us+5a ■ ■ ■ Us+L — 25aUs+L — 6aGs 

W[Cs+Sa] =UlW[Cs\Us. (A.3) 

But Us is just a unitary matrix, which means that if 9s+5a diagonalises W[Cs+Sa-] then Us9s+S(t 
diagonalises W[C's] with the same eigenvalues. Thus, with 9s = Ughg+Ss, 

9lW[Cs\9s = 9l^,^W[Cs+Sa]9s+5a = (A.4) 

The result follows by induction. □ 

Lemma Appendix A.3. Suppose that there is a field 9^ which diagonalises TT[C'cr] at each 
location a along the curve, with non-degenerate eigenvalues given by the diagonal elements of 
then 9\Ua9cr+5a is also diagonal. 

Proof We have, using the results of the previous lemma, 

=9l^,^W[Cs+Sa]9.s+Sa 

=9l^sM(^s9lW[Cs]9s9lUs9s+sa 

=9l^sMesU^^<^-^'9lUs9s+s. (A.5) 
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which means that 


= 0, (A.6) 

which, unless hF[C's] has degenerate eigenvalues, is only possible if 6 \Us 0 s+ 5 a is diagonal, as 
required by the lemma. □ 

There always exists a 9 field (unique baring the 17(1)^'^ phase factors and ordering of the eigenvec¬ 
tors) which diagonalises a Vb[Cs] G SU{Nc) as long as W does not have degenerate eigenvalues. 
From lemmas [Appendix A.l| and [Appendix A.3[ we see that this field is also a unique (up to 
the same caveats) solution for a field that diagonalises U. Thus we have shown that there is a 
SU{Nc )field 9 which diagonalises the gauge links along the closed Wilson Line. 

Appendix A.2. The Abelian Deeomposition in the eontinuum 

Here we construct an explicit form for A. The argument here is based on that presented in 




Lemma Appendix A.4. For any Field Y in the adjoint representation of a SU(Nc) He algebra 

^ = E [^^^]]) . (A-7) 

j ^ ' 

where are the diagonal elements of the lie algebra normalised so that tr(A'^ A^) = 26^^. 

Proof Y can be decomposed as H = A“T“, where A“/2 are the generators of the gauge group. 
A“ are Hermitian, traceless matrices which satisfy trA“A^ = 25°'^. We adopt the convention 
that the diagonal Gell-Mann matrices are referred to with the index while A“,A^,... 

may refer to any of the Gell-Mann matrices (including the diagonal components). The diagonal 
generators commute with each other, [A-^ , A^] = 0, and can be parametrised as diag(l, —1, 0,0,...), 
-^diag(l, 1, —2,0,...) ... —p=^==diag(l, 1,1,1,..., —(Ac — 1)). The other elements can be 

expressed in terms of the basis i?±Q-, where E^\a.\ consists of one element of 1 above the diagonal 
and its symmetric counterpart below it, while E_\a.\ has —i above the diagonal and i below it. For 
example. 


El =X^ = 

/ 0 1 \ 
1 0 ... 

E_i =X^ = 

/ 0 ••• \ 
i 0 ••• 


(; ; ■..) 


U : •••; 


Now, as is diagonal, [A-^, Ea\ must have the same non-zero components as Ea, and thus be a 
linear combination of Ea and E-a- Given that tiEalX^, Ea] = 0, we can write [X^,Ea] = af^E-a 
for some coefficient a which needs to be determined (a could be 0). Similarly, [Ea,E-a] is diagonal 
and traceless, so [Ea, E-a] = J2k some coefficients b. The following results then proceed 

immediately using the trace theorem ir{A[B, C]) = tr(C'[A, B])-. 

[A^[A^^;„]] =aialaEa 

tr(A„[A^ [A^ F;„]]) =2aiala = ti{[X^, Ea][E a, X^]) = -2iaif 

np = — np 

Uq, a 

6 ^ =\ir{X>^[Ea,E_a\) = ]Ar{E_a[X\Ea\) = ai 
tTi[Ea,E.ar) = - 8 = tviX^VaX’^bi) = 2^(6^)^ = 2^(a^)^. (A.9) 

k k 
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The first equality of the last equation follows from the specific form of Ea and E-a'- [Ea,E-a] 
has two non-zero terms on the diagonal, one 2i and the other —2i. We now have the necessary 
ingredients to prove equation (IA.7I) . Expanding Y = -\- daEa^ the right hand side reads 

^ QAnr(A^(c'=A'= + d„£;„)) + i[A^[A^c'=A'= + d„£;„]]^ (\i-l{ai^fd^E^ 

3 ^ j ^ ' 

=(d\^ Y daEo, = Y (A.IO) 

□ 

Remark In the Lemma above, we used a natural representation of and E^ in terms of the 
Gell-Mann matrices. However, it is easy to see that the same commutation relations between 
the generators apply if we rotate the basis X^ —^ 9X^9\ Ea —t 9Ea9\ with 9 G SU{Nc), 
while the traces of powers of the operators are also untouched. As the proof of (IA.7I1 just 
depends on the commutation relations and the traces, it is clear that a form of this equation 
is also valid in the rotated basis. In particular, this means that we can replace X^ in equation 
(IA.7I) with the rotated basis n\ leaving us with 

r = ^ QnHr(nJy) -F . (A.II) 

3 ^ ^ 

Using equation (lA.llll . we can use the defining equations to construct A. From the second equation 
da, we immediately have 

tr(n'’A) = tr(n'^A). (A.12) 

From equation (fT^ . 

0 =d^,v? - ig[A^,nA 

0 = ^ {[n\d^nA - ig[nA 

3 

= '^[nAdfj_nA + 2ig | 2A^ - ^ n-^tm-^A^ 

3 \ 3 

A^=Y1 Af,) + ■^[n\df,nA'^ . (A.13) 

The second field may be constructed from = A^ — A^. 

Theorem Appendix A. 5. The definitions of the restricted gauge field given in equations fa 
and i22\} . i.e. 

gAfj, = - + i9dfj,9^ 

Af,=Y^ QnHr(n^A^) -h 

are equivalent, where uf is a function of 9 and which needs to he determined. 

Proof We start from equation (IA.15I) . and write, with the help of equation (lA.llI) and = 


(A.14) 

(A.15) 
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= ^nHm^A^ + ^ [n^,0dfj,9^] 

= ^nHT{n^A^ - -68^6^)) + -68^6^ 
^ 9 9 


and if 

K - - ie8^e'<)), 

we recover equation (|A.14I) . 


(A.16) 


(A.17) 

□ 


Theorem Appendix A.6. The restricted field strength ean be written in the form given in 
equation \2t^) 

= \ {n^d^ti{n^ Ay) - 8yti{n^ A^)'^ + ^n'^tr(n'^ dj.n'']) (A.18) 

Proof We first of all construct F^y. For any field X, we have 

D^[A]Dy[A]X = 8^8yX - ig8f,[Ay,X] - ig[A^,8yX] - g^A^, [i,, A]] (A. 19) 

Anti-symmetrising this in the indices /i and v gives 

[D^[A],Dy[A]]X = - ig[8^Ay - 8yA^,X] - g\[A^, [A,, A]] - [A,, [A^, A]]) 

= - ig[8^Ay - dyAf, - A^], A], (A.20) 

and 

A^^[A] = 8f,Ay - 8yA^ - ig[Af,, Ay] (A.21) 

where we have made use of the Jacobi Identity [A, [B, C]] + [B, [C, A]] + [C, [A, B]] = 0 and we 
have used the definition of F^^y given in equation (fT^ . 

We may now use equation (IAA3I) to obtain an expression for F^y [A]: 

A^^[A] =i5^(n'^tr(n'^Ai,)) - ^8y{nHr{n^Af,))+ 

diJ.^[nAdyn^] - dy^[n^,8f_tn^] + 

i[nAr(n'’A^), [n'^,8yn’‘]] - ^ [n'’tr(n'^ A^,), [n'", d^n'']] + 

■^[[n\8f,n^, [n’^,8yn^]] 

^^riA8^j,tTn^Ay - dytTn^Afj,)+ 

{n^Ay){8 

itr(n'^A^)(a^n'^ - , In'",8yn’‘]])+ 

-^{8[8f_tn^ ,8yn^] + [[n^,8^^n^], In'", 8yn'']]). (A.22) 
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For the term proportional to trn^A,^, we can use the Jacobi identity and equation (lA.llI) to write 


[n\d^n^]] =d^n^ + ^[n\ [d^n^n^]] 

=d^n^ - ^[n^, [n^,df,n^]] 

=d^n^ — + in^tr(n^9^n'^) = in^tr(n^[n'^, 09^0^]) = 0 (A.23) 

For the term with commutators between Sj/U and d^n, we use d^n = [n, 09^0^] and expand 

S[d^n\d^n^]+[[n\dfj,n^, [n'',9^n'']] = 8[9^n^9^n^] + [[n\ [n-^,09^0'*']], [n^, [n'^,09i.0'^]]] 
=4(9^[n'’,9j.n'^] - 9i,[n'^, 9^71-^])+ 

4[n'’tr (71-^09^0^) — 20d^9\n^tr(n^6dv0^) — 209,^0^] 

=4(9;i[7i'^, [71-^,091,0^]] - 9i,[77-^, [n\9df^9^])+ 

4[77^tr(7i^09^0^) - 209^0^r^'=tr(7l'=09^0^) - 209^0^] 

= — 8[di_i{nHT{n^9d^9^)) — d^{nHr{n^ 9dfj,9^)) + 2dy9dfi9^ — 2dfi9du9^) — 
8(tr(77'’09^0^)9i,77'’ — tr(n-’09i,0^)9^77'^ — 2[09^0^, 09i,0^]) 

= — 8n^tr(—91/71-^09^0^ + 98,^9^ — 8^98^9^ + 8fi98u9^) 

= — 8n-^tr([n-’, 09^0^]09j,0^ — [n\98u9^]98^9^ + 77-’[09,y0^,09^0’^]) 

= — Sn^tr ( 77 -^ [09^0^,09^0^]). (A.24) 

This gives 

F^v[A] = ■i7i-^(9^tr7i-^A,, - 9^trn-^A^) - ■;^7i-^tr(77-^[09^0'l’,09^0^]). (A.25) 

z Zg 

Finally, we use tr(77^[77-^,X]) = 0 for any X, which follows from [ 77 ^, 71 -!] = 0, to show that 

tr77^[9^r7-^ , 9^77-^] =tr(r7^[[77-^ , 09^0^], [ 77 -^ , 09i,0^]] 

= - tr(77'"[r7-^, [09j,0'l', [ 77 -^, 09^0^]]]) - tr(77''[09j,0^ [[ 77 -^, 09^0'^'], 77 -^]]) 

= — 2tr(77^[09j,0^, 77^tr(77-^09^0^)]) + 4tr(77^[09i,0^, 09^0^]) 
=4tr(77''[09^0t,09^0t]), (A.26) 

which gives the alternative form for the restricted field strength tensor 

Ffj,^[A] = \n^{8f^tvn^Au - 8Am^AfAj + ^77-^tr(77^ [ 9 ^ 77 '“, 9„77'']), (A.27) 

2 8g 

in agreement with equation (IA.18|) . The factors of 1/2 and 1/8 depend on our normalisation 
convention for the generators of the gauge group. □ 

Appendix A.3. Exponents of the Lie algebra 

In this section, we give a few results which were needed in section [5l 

Lemma Appendix A.7. For non-commuting objects X and Y in some Lie algebra, 

e^Ye-^ = r + [A, r] + i [A, [A, F]] + ... (A.28) 

Proof The result is standard; it follows by expanding in a Taylor series and commuting A 
term by term through Y. □ 
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Lemma Appendix A.8. For non-commuting X and Y in some Lie algebra, 


(.X+Y ^ gXgYg-i[X.F]gi([X.[X.F]]-2[Y.[X.F]])g-i([X,[X,[X,Y]]]+3[[[X,y].X,Y]]+3[[[X,F].F].F]) 


Proof See, for example, |4Cll |. 


(A.29) 

□ 


Lemma Appendix A.9. For non-commuting objects X and Y in some Lie algebra, where X 
is differentiable with respect to a variable s and we can neglect terms of OfSs'^) 


^X{s) SsY -Xis-eSs) ^ 


(A.30) 


where we use the notation 


[A-F]^[A,[A,...[A,r]]]. 


(A.31) 


n terms 

Proof Using a Taylor expansion, and the results of equations (IA.28|) and (jA.29l) . we have 

^Xis)^SsY^-Xis+Ss) ^ ^ ^^Y^^_xis)-5sO,Xis) ^ 

= (1 + SsY + J2 ... + OiSs"^) 


n>0 




•>0 ^[^"'5sY]-5sa.x-E„>i o{6s'^) 


where X = X{s). Equation (IA.30I) follows immediately. 
Lemma Appendix A.10. For <p, X G su(2). 


(A.32) 

□ 


X . E (X - (A.33) 

Proof We may always parametrise (j) = for y S SU{2)/U{1) and real a. This means that 

equation (lA.llI) is applicable, and we have 


[if), [i(^. A]] = 2(^tr((/)A) — 2tr(j)‘^X. 


(A.34) 


Therefore, for integer n > 0, 

[(#)2"A] =(-2tr(/)2)”A + {-2tT(l)^)^-^2(j)tT{(j)X) 

[{i(j)f^+^X] =i{-2iY(jf^)^[(j),X]. (A.35) 


Therefore 


X + [i(j), X\ + ^ 


n>0 


L(2n)! 


[{icfY^x] + 


1 


(2n+ 1)! 




cos \/ 2iY(jF' i X — cj) 


^tr((5iA)^ ^ .sin-\/2tr^ 

tr02 J * y^tr^ 


[cf,X] + cf 


ti{(j)X) 
trcjp ’ 


(A.36) 


which is the result we wanted to prove. 


□ 
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Appendix B. Area Law and the String Tension 

We expect that the average number of the merinthons contained within a loop will be propor¬ 
tional to the area of the loop. Given that different configurations have statistically independent 
numbers of merinthons, we may model the distribution of the number, n^, of merinthons of a type 
j (i.e. merinthons constructed from A-^), within the Wilson Loop across the configurations using 
a Poisson distribution with a mean value /ij which is the density of merinthons multiplied by the 
loop area. Each merinthon (or merinthon -anti-merinthon pair) will contribute a certain value 
to the Wilson Loop, and this will be distributed according to some distribution P{x^) which 
is an even function of ^is positive winding number and negative winding number merinthons are 
equally likely. Neglecting any perimeter contribution from non-merinthon effects, the expectation 
value of the Wilson line will be 

tr((0ttE0))=trJ]^^e-'^^ iyj dx’^P{x\x\ ■ ■ 

3 ‘^3 ^ 

=tr n (g 

3 

We may expand the exponential in terms of trigonometric functions of x3 ■ For example, in 

SU(2), 

X = cosy^ + Jsiny^A^, (B.2) 


and in SU(3), 


s*"V+aV = if2cos 


2y« y« 2y« 

—= cos y"^ -I- cos —-= d 2i sin cos y"^ — f sin —^ 
Vi V3 ^3 \/3 


X 


-k A-' ( -sin-^smy 


3 I • X • 3 

' I COS —= Sin Y 


f X® 3 2y® . . X® 3 ■ • 2y 

( cos cos y — cos -I- i sm cos x sm 


Vi\ Vi 


Vi 


Vi 


(B.3) 


In each case, the imaginary terms are eliminated by the integral over y, so as long as the distiribu- 
tion P is not too strongly peaked at y = 0, we are left with a result which resembles ^ 

where ^ is proportional to the area of the loop (and depends on the number of structures) while a, 
which represents the quantity obtained after the integral over y, is independent of the area (and 
depends a single structure); each diagonal term exponentially decreases with the area of the loop 
and the trace will be dominated by whichever term decreases the slowest. This gives an area law 
scaling for the Wilson Loop. 


Appendix C. SU(3) Merinthon algebra 


The following commutation relations are useful in evaluating the product of 9 matrices in 
SU(3). ^ 1 , 4 )2, ... are defined in equation 


[4>i, =2i4>i 

[,(,l,^/3A®] =0 

[4>i,42] = - i cos(ci -I- C 3 - C2)^3 “ * sin(ci -I- C 3 - C2)4>Z 
[4>i,43] = - i cos(ci -I- C 3 - C 2)^2 + * sin(ci -I- C 3 - C2)4>2 
[<(>1) ^i] = ~ 2iA® 

[</>l, 02] =f COs(ci -I- C 3 - C2)03 “ Z Sin(ci -I- C 3 - C2)03 


[01, [01, A®]] =4A® 
[0i,[0i,v/3A®]] =0 
[ 01 , [ 01 , 02 ]] =02 
[01, [01,03]] =03 
[01, [01,0l]] =401 
[ 01 , [ 01 , 02 ]] =02 
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[01,03] =iCOs(oi -1- 03 - 02)02 -1- isin(oi -|- O 3 - O 2)02 

[01, [01,03]] =03 


[02, A®] =i02 

[02, [02, A®]] =2^ 

® + A3 

2 

[02,y3A®] =3i02 

[02,[02,V3A®]] = 6 ^ 

® + A3 

2 

[ 02 , 0 i] =* cos(oi -1- 03 - 02)03 + * sin(oi -1- O 3 - 02)03 

[ 02 , [ 02 , 0 l]] =01 


[02,03] = - * cos (03 -1- 0 i - 02)01 “ * sin(oi -1- 03 - 02)01 

[02, [02,03]] =03 


[ 02 ,0i] =i cos(oi -1- 03 - 02)03 -1- i sin (02 - oi - 03)03 

[ 02 , [ 02 , 0 l]] =01 



[02, [02,02]] =402 


[ 02 , 03] = - * cos(ci -1- 03 - 02)01 -1- * sin(oi -1- 03 - 02)01 

[02, [02,03]] =03 


[03, A®] = - 103 

[03, [03, A®]] =2^ 

\/3A® 

2 

[03,V^A®] =3i03 

[03,[03,y3A®]]=-6^' 

[ 03 ,0i] =*cos(oi -1- O 3 - 02)02 - zsin(oi -1- 03 - 02)02 

[03, [03,0l]] =01 


[03,02] =i cos(oi -1- 03 - 02)01 -1- i sin(oi -I- 03 - O 2 ) 0 i 

[03, [03,02]] =02 


[03, 0l] = - * COs (03 -1- Oi - 02)02 “ * Sin (03 -1- Oi - 02)02 

[03, [03,0l]] =01 


[03, 02 ] = - * COs(oi -1- 03 - 02)01 -1- i Sin(oi -I- 03 - O 2)01 

[03, [03, 02]] =02 


r .1 7 1 n A3-\/3A® 

[03,03] =2j 2 

[03, [03,03]] =403- 

(C.l) 

We obtain, using eouation (IA.28L the definition of 0 given in 

equation (1431) and the commutation 


relations listed in equation ( 1 C. 111 . 


_gj</i3a3gi02a2 gj^ioi ^3g-i0iai ^-ic/>2a2 g-i03O3 

cos(2ai) — 01 sin(2ai) 

sin( 2 a 2 ) 7 cos( 2 a 2 ) — 1 


_gj03a3gi02a2 


g —10202 g —*0303 


0*0303 


cos(2ai) — 


-02 -I- 


V3A' + + I - 


2 ^ 4 4 

sin(2ai)(cos(a2)0i — sin 02(03 cos(c2 — ci — C3) -I- sin(c2 — ci — C 3 ) 03 )) 

= — cos( 2 ai)^-^-^ sin 03(01 cos(c2 — C3 — ci) — 0 i sin(ci + C3 — C2))— 


0 — 10303 


,sin 2 o 2 7 
C 0 S( 2 oi)- - -cos 0302 - 1 - 


cos(2ai)- 


cos2a2 —1 f sin(2a3)- 


- 3 - 


-03 + ^A® 


1 -I -3 003(203) 3 


-l- — A ®(1 — 003203) ] -|- 


cos(2ai) 


cos 202 -f 3 3 003(203)-I- 3 ^ 1 — 003(203) , sin(2o3) 7 


A® 


03 - 


4 4 2 

sin( 2 oi) cos02(02 sino3 sin(c3 -I- ci — C2) -I- 02 sinos 003(03 -l- Ci — 02) -|- 0 i 00303)4- 
sin(2oi) sino2 cos(c2 — ci — 03)03-!- 

sin(2oi) sino2 sin(o2 — ci — 03) ^cos(2o3)03 — — \/ 3 A®)^ . (C.2) 


We also obtain, 

VSn® =e“303g*O202y3;^8g-*a202g-i03a3 
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_gia3<^3 


A' 


3cos(2a2) —3 /-,o3cos(2a2) + 1 3111 ( 202 )- 

'-^-+ v'3A -^-3--- <j )2 


3 - 3 - 

= - - siii(2a2) sin(o3)(^i cos(c2 - C3 - ci) - sin(ci -|- C3 - C2)) - - sin(2a2)^2 003(03)-!- 


3co3(2o2) — 3^^3003(203)-I -3 ^ 1 — 003(203) ^ sin(2o3)- 


-(A^ 


h)+ 


+ + 73(1 - cc,s( 2 a 3 ))). 


A 3imilar oaloulation give3 


o^d^e =Si^ + 82 ^ + ^3M + 


(C. 3 ) 


(C. 4 ) 


where 


_^—idiX^ — idsX^Q gidsA^+idgA® 

=iX^d^d 3 + iX^df,ds, 


(C. 5 ) 


_^-id 3 X^ -idsX^ q ^^iai(f>i ^^id^X^ -\-idsX^ 

i(j)id^ai + i9^ci(cosai sinai^i — sin^ aiA^) 


_g— —id^X^ 


pZdsA^+idgA® 


(C.6) 


_^-id3X^-idsX^ ^-iai(pi^-ia2(p2Q f^^ia2<p2^^iai(f)i ^id3X^-\-idsX^ 


_^—id^ A —ids A® 2(^1 ai 


dfia 2 i 4>2 + 19^02(003(02) sin(o 2)(()2 - (A^ + a/SA®)) 


^i(piai gidsA^+idg A® 


_^—id^ X^ — ids X^ 


2 

— »9^C2 (a/SA® — oo3(2oi)A® — 3iii(2oi)^i) + 


id^a 2 {cosai(j )2 - 3in01(003(02 - oi - 03)^3 - 3in(o2 - ci - C3)(/)3))+ 
*9^02^^^^^^(0030102 + 3inoi(003(02 - 0i - 03)93 + sin(02 - ci - 03)^3)) 


pidsA'^+zdg A® 


(C. 7 ) 


and 


g-io2 02 e“*“303^^ (^gia3 03 ^g*O2 02 giaiipi gidsA^'+idg A® 


_g—ids A —ids A g—iai^i g—ia2 02 


gio2 4>2 giai 01 gids A^ +id8 A® 
_g—idsA® —idsA® g—iai0i 


d^a3i(j)3 + id^cs foo 3 033100393 - (A^ - v^A®' 


*9^03 (93 003(02) - sln(a2)(9i 003(02 - 03 - oi) + 9i sln(oi + 03 - 02))) + 

idfj,C3 COSO3 slno3 (93 003 02 — 3ln02(91 cos(oi + 03 — 02) — 9i sln(oi + 03 — 02)) — 

• 2 

*9^03_ 1/3A® - (003(202) - 1 )(a/ 3 A® + A®) - 310(202)92) 

giai 01 gids A®+id8 A® 


_g—ids A®—ids A® 
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Appendix D. Numerical Methods 

Appendix D.l. Algorithm for fixing 9 

It is simplest to find 6 by solving for the eigenvectors of the Wilson Loop; however we employed 
a different technique (the method discussed below has the advantage that it can be easily adapted 
to other rules for choosing 9 discussed in the literature than just the one we use in this paper and 
to other problems such as gauge fixing and solution of the defining equations. It is also reasonably 
quick and efficient.). 

We wish to find the 9 that solves 

= (D.l) 

for each diagonal and gauge link along the Wilson Loops. This is equivalent to solving 

0 = _tr[([/^,^n^_^^ - niUx,p){Ul^^ni - n^+^C/],^)], (D.2) 

which is solved by the 9 which minimises 

F[9] = E [ 4 ^ - ’ (D.3) 

j 

with 

ni = (D.4) 

and Al = A^ or A® and 9 G SU{3) (in fact, the solution we want is at F[9] = 0). The trace is over 
both lattice sites and SU(3) indices and V is the lattice volume. 

We solve this minimisation problem by combining two different algorithms; one based on 
molecular dynamics, which is an extension of the steepest descent method, and the second a 
Newton-Raphson algorithm. The goal is to have the molecular dynamics solve the system to a 
sufficient accuracy that the Newton-Raphson algorithm will converge, and then Newton-Raphson 
rapidly finds the solution to a high precision. The combination of these two methods, while it 
may not be the most efficient algorithm, worked well on all the lattice volumes we used. We 
also employed similar algorithms to gauge fix and solve the defining equations of the Abelian 
decomposition. 

Having found 9, one must subsequently order the columns of 9 and fix the U(l) phases according 
to the chosen fixing condition. We ordered the columns according to decreasing real part of 
eigenvalue, and either fixed 9 using c? = 0 or by having the same value for each u around the 
Wilson Loop, depending on the observable. 
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Appendix D.1.1. Molecular dynamics 

For the molecular dynamics algorithm, we introduce a momentum field tt{x) conjugate to 0 
and define a fictitious energy 

E = ^tiTT^ + F[e]. (D.5) 

TT here represents a traceless Hermitian matrix on each lattice site. 

We evolve d in a fictitious computer time r using 

-^Ox = iOxT^x, (D.6) 

ar 

or, as 6t —>■ 0, 

Oxir + 6 t) = 9x{r)e^^^^^. (D.7) 


The second equation of motion is taken from the conservation of energy. 


d d 

0 =tr7r—TT + —F 
dr dr 

F =Retr(7r(a;)G(a;)) 


dr 

G' =i 


olUx.udx+uX^oi^Ml jxX^+elul_, jx-uX^0L,ul_ 


.9xX^ 


X^0lul_^Jx-pX^9l_-Ux-y,^9x + X^9lUx,,9x+pM9l^^Uljx 


(D.8) 


G is taken to be the hermitian, traceless part of G', and in our particular problem (since the 
and A® components of 9 are arbitrary) we also set the diagonal components of G to be zero. We 
then evolve 9 and tt using the Omelyan integrator [^, for a time step St 


'^x ^'^x (ySxGx 

9x ^9xF^^- 
T^x -tTTa; - (1 - 2a)5TGx 
9x -^9x0^ 2 

TTx ^Tx - aSrGx- (D.9) 


This conserves energy up to 0(<5t^). a should be tuned to have the best conservation of energy 
for a given dr, while St should be the maximum value where the algorithm remains stable (i.e. 
energy is conserved). We repeat this integrator a number of times to form a trajectory, ending the 
trajectory when F[9] stops decreasing. We select the 9 which had the lowest F[9] over the course 
of the trajectory, and repeat until F[9] is below a target threshold. 

The trick is to start each trajectory with tt = 0. This means, at least for the first few steps, 
ti'TT^ must increase, and since energy is conserved, F must decrease. While this method cannot be 
used to find a precise solution for the minimisation of F"[0], as the force G becomes zero close to 
the minimum and the method slows down, it rapidly converges towards a solution, and gets close 
enough to allow the use of more powerful methods. The initial step is equivalent to the steepest 
descent method; but as the trajectory progresses we incorporate information from the forces at 
larger distances from the lattice site. 

The exponential of the anti-hermitian 3x3 matrix is calculated using the Caley-Hamilton 
theorem, following 341. 


Appendix D.1.2. Newton-Raphson 

Once we have found the minimum to a sufficient accuracy, we employ a Newton Raphson 
algorithm to polish it to an arbitrary accuracy. 
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We define a force according to 


ra_dmn .a 


de, 

= — Re tri 
— Re trz 




(D.IO) 


This force is related to the molecular dynamics force above by G® = itr(A“Ga;), and is essentially 
a vector representation of the force rather than a matrix representation. The force gradient is 
given by differentiating F twice 


Hx,y = - 4.s/Re trz 


— <5£!;,j/Re trz 


- 6x-fi,,yRe trz 


- (5a;+/i,yRe trz 

We then update 6 according to 


Ux-Mx,, [A“, x^]]9lul_^jx-yX^el_f, 
Ux,y9x+f.X^9l^^Ul^9x[Xh, [A“, X^]9l 
Ux-MX% X^]9lul_^Jx-y[X\ X^]9l_^ 
Ux,y9x+f,[X\X^]9l^.Uljx[X^,X^]0l 


-'X ' ^x' 




(D.ll) 


(D.12) 


In this particular problem, the sum over a and b excludes the indices related to the diagonal Cell 
Mann matrices (where G is zero, and this removes the zero eigenvalues from H). H^'y can be 
easily inverted using a CGNE algorithm with odd/even preconditioning. 

This is equivalent to the usual Newton-Raphson minimisation algorithm, and has the usual 
advantages and disadvantages associated with Newton-Raphson. It has superlinear convergence 
as long as H^'y contains no small eigenvalues and we start close enough to the desired solution. 
Restricting the sum over A“ to the off-diagonal elements ensures that the inversion of H^^y is 
smooth enough, and we did not encounter any difficulties: F[9] was reduced to zero within machine 
precision in two or three iterations. 

Once we have found a 9 for which f = 0, it is straightforward to fix the ordering of the 
eigenvalues of the Wilson Loop and the U (IVc— 1) phase according to the choice of hxing condition. 
We sorted the eigenvectors according to their real part, and fixed the phase so that 9^U9 is equal 
for all the links along the Wilson Line (so under a gauge transformation 9 ^ A9 without any need 
for any additional gauge fixing). 


Appendix D.2. Numerical solution of the defining equations 

Given zzj, we wish to find the solution to the lattice defining equations 

u jA — T7^ —n 

ztr(zz^(X^-(l^)t)) =0, (D.13) 

for Ux,y = Ux,fiXly. The solution of these equations is given by the X € SU {Nq) which minimises 
the functional 


g[Xy] =gAXy]+gAXy] 

92[XA =E [-AAni{Xx,y - XlyMtriniiXx,^ - n,y))) 

3 


3 


{Ux-y,yXlyniXx,yUl_y^ 




— n 


3 

X — fJ. 


) 


2 


(D.14) 
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for each direction /i. Each can be solved independently, although in practice we set up our 
algorithm to solve them all simultaneously. Clearly, g[X] > 0, with the equality achieved when 
the defining equations are satisfied. If there is an exact solution to the defining equations, this 
procedure will find it. If not, we at least achieve something close to the lattice solution. In practice, 
by varying the initial guess for X we found that there were usually three or four solutions for each 
Xfj,. We choose the solution with the largest trX, a condition which is both gauge invariant and 
satisfied on those links along the Wilson Loops used to define 9, where U = U and X = 1. 

The same minimisation routine is used as in the previous section (solving for X), combining 
molecular dynamics with Newton-Raphson minimisation, and while we make no claims that this 
is the most efficient algorithm, it proved to be efficient enough on our configurations. There were 
no issues concerning algorithmic stability or convergence. 

The Newton-Raphson force and force Gradient are 


= - 2*(tr«(l,,^ - Xl^))){t,{nUX^,^X'^ + A“Xt_^)))+ 


2 ztr 


„Xl „niX^ „Ul 




(jX^x-fi,fiXl^^ni.Xx,nX Ux-fi,fiX Xl^^ni.Xx,fiUl_^^^ 


=2(tr«(l,,^A*' + X'^Xl^mtv{ni{Xx,^X^ + A“lt_^)))+ 
2(tr(ni(l,,^ - Xl^))){t^{ni{Xx,,X'>X^ - X^X'^Xl^)))- 


2 tr 


2 tr 


2 tr 


t 


(Ux-fi,fiXl^^ni.Xx,fiX — Ux-ij,,ij.X Xl^^ni.Xx,fj.U, 

(Ux-fi,fiXl^^ni.Xx,fj.X — Ux-fi,fiX Xl^^ni.Xx,ij,Ul._^^^ 

(jJx-fj,,fj,Xl.^^ni.Xx,fj,Uj._ij^i^ 

(^~Ux-fj.,fj.X Xl^^ni.Xx,ij,X + Ux-fi^fiX X Xl^^ni.Xx^fj,ul_^^^ 

(jJx-fj,,fj,Xl.^^ni.Xx,fj,U}._ij^ij^ — n: 


3 

X — fl 


(Ux-ii,iiXl.^^ni.Xx^fj,X X Ux-fj,,fj,X Xl.^^ni.Xx^fiX 


(D.I5) 


and the molecular dynamics force can be constructed from the Newton-Raphson force using Gx = 
A“G“. 
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